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ABSTRACT 

A COMPUTER MODULE USED TO CALCULATE THE HORIZONTAL 
CONTROL SURFACE SIZE OF A CONCEPTUAL AIRCRAFT DESIGN 

Stephen Mark Swanson 
January 1990 

This paper discusses the creation of a computer module used to calculate the size ot 
the horizontal control surfaces of a conceptual aircraft design. The control surface size is 
determined by first calculating the size needed to rotate the aircraft during takeoff, and, 
second, by determining if the calculated size is large enough to maintain stability of the 
aircraft throughout any specified mission. The tail size needed to rotate during takeotf is 
calculated from a summation of forces about the main landing gear ot the aircraft. 1 he 
stability of the aircraft is determined from a summation of forces about the center of 
gravity during different phases of the aircraft's flight. Included in the horizontal control 
surface analysis are: downwash effects on an aft tail, upwash effects on a forward canard, 
and effects due to flight in close proximity to the ground. 

Comparisons of production aircraft with numerical models show good accuracy tor 
control surface sizing. A modified canard design verified the accuracy of the module tor 
canard configurations. 

Added to this stability and control module is a subroutine that determines one of three 
design variables for a stable vectored thrust aircraft. These include forward thrust nozzle 
position, aft thrust nozzle angle, and forward thrust split. 
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CHAPTER 1 


Introduction 

In the aircraft conceptual design process, there are five major areas in which the 
designer allots most of his or her time and effort. These include aircraft layout, 
aerodynamics, weights, propulsion, and performance. In each of these areas the designer 
goes through a design process that includes a large and complex senes of decisions and 
calculations to determine the design parameters of the aircraft. After initial parameters 
have been determined, the design is compared to any specified requirements, appropriate 
changes are made, and then another series of decisions and calculations is completed to 
refine the design. This cycle is repeated until the aircraft design created meets the specified 

requirements. 

Once the conceptual design has become refined, the designer turns to more detailed 
areas in the aircraft design process. One of these areas is that ot stability and control, 
where the sizes of the control surfaces on the aircraft are determined, the stability of the 
aircraft is determined, and where a more refined estimate of aircraft components weight 
and position is completed. Once again, as in the conceptual design areas, the detailed 
design areas involve a cycle of calculations and decisions in order to create a more refined 

aircraft design. 

NASA Ames Research Center has created a computer program which does the 
calculation part of the conceptual design series. This design program, called ACS Y NT 
(AirCraft SYNThesis design program), helps the designer in the five areas ot the 
conceptual design. This includes geometric layout, aerodynamic analysis, weight 
calculations, propulsion estimates, and performance analysis. This program allows the 
designer to concentrate more on the decisions and less on the complex calculations. 
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To enhance the conceptual design process of the ACSYNT design program, the 
decision was made to create a module to perform calculations necessary for a longitudinal 
stability and control analysis. By performing more analysis of the aircraft's design at an 
earlier stage, a more refined product is created during the conceptual design, reducing the 
amount of effort required to finalize the aircraft design. 

This paper discusses the development of the stability and control module. This 
module will be used in conjunction with the ACSYNT design program to enhance the 
conceptual design process. The module is used to determine the aircraft's center ot gravity 
by positioning the different components of the aircraft at positions specified by the 
designer. It is used to determine the shift in center of gravity due to fuel and weapons 
usage during the aircraft’s mission. Finally, it is used to calculate the horizontal control 
surface size needed to maintain controllability during take off, mission completion, and 
landing. The module is used to determine the horizontal control surface size ot an aircraft 
which uses either a conventional aft tail, a forward canard, or both. 

In order to determine the stability of the aircraft, there are several aerodynamic 
parameters that the module must determine. These include the shift in the aerodynamic 
center of the wing with respect to Mach number, and the lift curve slope ot the canard if 
one is being used (lift curve slopes for the wing and aft mounted tail are calculated using 
ACSYNT's aerodynamics module). Also, the downwash of the wing on the tail, the 
upwash of the wing on the canard, and ground effects are determined using the module. 

Also included in this module is an analysis that uses the stability conditions during 
landing to calculate three vectored thrust design parameters. The analysis includes the 
position of the forward thrust vector, the angle of the aft thrust vector, and the amount ot 
thrust split between the forward thrust vector and the total thrust. This allows the designer 
to create a vertical landing aircraft design which is stable during transition from torward 


flight to hover. 



CHAPTER 2 
Center of Gravity 

One of the more important parameters used in stability calculations is the center ot 
gravity (C.G.). The C.G. is the equilibrium point where the weights of the different 
components act as one force at one point. It is a reference point in the aircraft s design in 
calculating both the horizontal control surface size, and in determining the stability ot the 
aircraft design. 

Determination of the C.G. is derived from the definition, finding the one point where 
all the weights act as one force. This is done by summing all the weights and moments of 
the aircraft components about the nose of the aircraft, and dividing the two as shown in 
Equation 1. 

Y Moments ( 1 ) 

C.G. = ^ — 

Weights 

The C.G. position is determined in units of length from the nose. The C.G. provides 
a reference point for the summation of aerodynamic moments of the aircraft. 

Since ACSYNT is a conceptual design program, there is little information calculated 
in ACSYNT about the positions of the different aircraft components. The program 
calculates only the positions of the engine(s), wing, vertical tail, aft tail, and forward 
canard. Positions for all other components are selected by the designer and input into the 
stability program. 

The component positions that ACSYNT does calculate are determined in the 
following manner. The position of the engine is determined from the value ot the length 
and position of the engine pod supplied through ACSYNT’s geometry module. This can 
be modified in the stability module by the designer through adding or subtracting a fraction 
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of the pod length. The locations of the wing, horizontal tail, vertical tail, and canard are 
determined from values of the quarter-chord (C/4) points for each aerodynamic surface as 
supplied by ACSYNT. The stability module allows the designer to vary these positions 
by adding or subtracting a fraction of the mean aerodynamic chord (MAC). 

All other component positions are left as the designers choice, allowing placement of 
the components at any point along the fuselage. The component positions are calculated 
using multiplying factors input by the designer. Equation 2 shows an example ot this 
process where the C.G. of the fuselage is calculated as a function of the fuselage length 

(BODL), and a multiplying factor (Xf fus ). Designer input multiplying factors allow for a 
more flexible design process, giving the designer more freedom in the internal layout ot 
the aircraft. If the designer does not desire to specify any or all ot the component 
positions, default positions inside the module are used. 

x cg| „ = Xf fus * BODL <2) 

When determining the stability of the aircraft, it is important to calculate the shift in 
C.G. throughout the night as fuel is used and/or weapons deployed. If the C.G. travels 
too far forward or aft during the flight, it can put the aircraft in an unstable condition, 
making it impossible to fly. 

The ACSYNT program allows the designer to divide the aircraft mission into a 
maximum of 12 different phases. The designer specifies for each phase the altitude, range, 
Mach number, and engine thrust setting of the aircraft. The designer can also specify it 
any weapons are to be deployed at the end of a phase. For each phase, ACSYNT 
determines the lift, drag, and fuel usage. It can also determine the optimum altitude, and 
optimum Mach number. ACSYNT then determines a new aircraft weight for the next 
phase of flight by subtracting out the weight of the fuel used and weights ot any weapons 
deployed. 
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To determine the shift in the C.G. position at the end of each phase, the weight sum 
and moment sum are calculated using the weights for each phase. To further enhance the 
capabilities of the stability module, modifications were made to include the weight ot any 
external tanks. When the amount of fuel used becomes greater then the amount of fuel 
stored in the external tanks, the module removes the external tank weights, removing them 
from future C.G. calculations. 

Once the C.G. had been determined and the C.G. range calculated, the position ot the 
C.G. with respect to the point of ground contact of the main gear is evaluated. The 
position of the point of ground contact is important when determining ot the horizontal 
control surface size necessary for take-off rotation. From U.S. Navy aircraft design 
specifications 1 , the C.G. should be located at least 15 degrees forward ot the main gear, 
as shown in Figure 1. This assures the aircraft will not accidentally tip over during 
ground handling operations, and that enough weight rests on the nose gear to maintain 
positive nose wheel steering of the aircraft. The stability module compares the C.G. 



1 Curry, Norman S.; Aircraft Landing Gear Design: Pri ncipl es and Practices , 
AIAA Education Series; American Institute of Aeronautics and Astronautics Inc. 
370 L’Enfant Promenade, S.W., Washington D.C., 20024; 1988; 43 - 50. 
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position calculated for each phase of flight to the position of the main gear. If the 15 
degree requirement is not satisfied, the module moves the point of contact aft a fraction ot 
the fuselage length (0.01 * BODL). It recalculates the C.G. until the requirement is 
satisfied. Comparison between the C.G. and the point ot contact during all phases ot 
flight assures the aircraft will meet the requirements if an inflight emergency were to 
require landing prior to mission completion. 



CHAPTER 3 

Horizontal Control Surface Sizing 


There are two different methods used in determining the size of the horizontal control 
surface. The first method determines the size of the surface necessary to rotate the aircraft 
during takeoff. The second method, discussed in a later section, determines if the control 
surface size calculated by the first method is large enough to satisfy the forward and aft 
C.G. limits of the aircraft. When both methods are satisfied, the aircraft will have a 
horizontal control surface large enough to maintain positive longitudinal control. 

The first method sizes the horizontal control surface by a summation of moments and 
forces about the point of ground contact of the main landing gear on the aircraft, as shown 
in Figure 2. The forces on the aircraft at takeoff rotation include: wing lift, horizontal 



Figure 2: Summation of Moments about the Main Landing Gear 


control surface lift, wing-body pitching moment, aircraft drag, aircraft weight, and engine 
thrust. Also included is the aircraft acceleration. The moment arms for each of these are 
taken with respect to their positions from the point of ground contact. This point is used 
since the aircraft needs to rotate about the point of contact during takeoff. The wing lift is 
positioned at the aerodynamic center (AC). It is modified to include ground effects, which 
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discussed in a later section. The aircraft's drag and acceleration are placed along the 
aircraft centerline (ACSYNT calculates the drag of the entire aircraft and uses that in all its 
calculations. This includes the drag of the wing, fuselage, and tail surfaces. The total 
drag force as defined by ACSYNT is therefore put at the aircraft's centerline). The 
maximum takeoff weight, as determined by ACSYNT’s weight module, is located at the 
C.G. The horizontal control surface lift, also modified for ground effect, is placed at its 
AC. In ACSYNT, the ACs of all lifting surfaces are the quarter-chord point of the mean 
aerodynamic chord (MAC). The wing-body pitching moment is positioned about the point 
of ground contact. For the takeoff analysis, a conventional takeoff is assumed. This 

means that the forward and aft thrust vectors (Tj and T 2 ) point in the alt direction. 
Therefore the only moments created by the thrust vectors are due to the vertical 

displacements (Z T1 and Zrpo) of the two vectors. 

The horizontal positions of the wing lift, aircraft weight, and horizontal control 
surface lift depend on whether the horizontal control surface used for pitch control is a 
conventional aft tail or a forward canard. For an aft tail configuration, the wing lift is 
placed forward of the C.G. and the point of ground contact, and the horizontal tail lilt is 
placed at its AC in the rear of the aircraft. For a canard configuration, the wing lilt is 
placed behind the C.G., and the canard lift is placed at its AC in the front of the aircraft. 
This leads to two different equations for the summation of moments about the point of 
contact of the aircraft. A third equation is derived for the case of both a canard and a 
horizontal tail. In this case, the canard is assumed to be used only for trim." Its size and 
shape are specified by the designer, who then uses the stability module to size the ait 


2 The ACSYNT program in its Aerodynamics and Geometry modules is unable to handle 
cases where there exists both a horizontal tail and a canard on the same aircratL Having 
the designer input the canard size and shape, allows the designer to evaluate those cas . 

where both exist. 
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horizontal tail. This uses the same equation as the case with a horizontal tail, but includes 
a term for the moment created by the lift of the canard acting at its AC. 

Examining first the case for which there exists a horizontal tail only, the equation tor 
the summation of moments created by the aerodynamic forces about the point ot ground 
contact 3 is seen in Equation 3. This is modified to account for any vertical thrust offset as 
specified by the designer. 


M mg = - C Mwb9 Sc - [ Z mg _Z T 1 ] T l ' [ Z mg- Z T^ T 2 ' [ X mg- X cg] W 

W . r 7 

+ [ X mg‘ X ac] L wb + [ Z mg' Z D] D + T VZ mg ' |^ X ac ht ' X mgi 


(3) 


L 


hi 


In this equation, the lift of the wing (L wb ) is replaced with the classic definition 


Lwb=C Lag a w qS 


(4) 


where the subscript g represents changes due to ground effect which is discussed in the 
section titled "Aerodynamic Calculations". 

From Newton's equation of motion F = ma, the aircraft acceleration component is 
replaced by the forces in the horizontal direction. These forces include thrust, drag, and 
the reaction force of the aircraft acting on the runway. This replaces the acceleration term 

for the aircraft with 


— V = T - D - uR 
g 


3 Roskam, Jan.; Airplane Flight Dynamics and Automat ic Right Controls, Part i; 
Published by the author; 519 Boulder, Lawrence KA; Second Printing, 1982; 373-3/5 
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From the summation of forces in the vertical direction, it is possible to solve for the 
reaction force of the aircraft acting on the ground. 

R = W - L w5 - L hl (6) 

In order to solve for the minimum tail size needed to rotate the aircraft, the forces 
about the point of ground contact must be balanced. To achieve this condition, M mg in 

Equation 3, is set equal to zero. The horizontal tail lift is solved using the 

combination of Equations 3, 4, 5, and 6, as seen in Equation 7. Equation 7 is used to 
calculate the necessary tail lift in terms of the different aerodynamic forces. 

L ] + L 2 - L 4 - L 6 + L 7 - L 8 (7) 


where 


Li - - C Mwb q S c 


U = 


X m g ' X C g + M- Zmg 


W 


u = 


Xmg' (X cg - Xj,) + M-Z mg 


L 7 — 


X mg " X aCwb +d Z mg 


-wb 




La = 


X'i’2"^' Xcg' X m g ■ ftZ m gj 


A = 


X a c hl ' X mg ’ 


It is important to note that the lift required by the horizontal tail to rotate about the point ot 
ground contact is negative in value (that is pointing down), which requires a negative 
horizontal tail angle of attack. From the classic equation for lift, the horizontal tail lift is 
expressed in Equation 8. 
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Lht = CL aht g (ttw + >ht - iw + T Se) 9 Sht - Clqiu e g q Sht ( 

Where the angle of attack of the tail as shown in Figure 3. Rearranging Equation 8, 



creates Equation 9. Equation 9 is used to solve tor the horizontal control surface size at 
take off rotation, while in ground effect. 

s Ut (9) 

C L ah t g k + ihi - >w + t 8 e ) q - C^, £ g q 

For the case where there exists both an aft horizontal tail and a forward canard, the 
assumption is made that the canard is used only for trim, and not tor longitudinal control 
purposes. With the size of the canard defined by the designer, the stability module is 
used to determine the size of the aft horizontal tail. Modifying Equation 7 to include the 
moment created by the canard lift placed at its AC, leads to Equation 10. 

L, + L 2 - L 4 - L 6 + L 7 - L k + L 9 (10) 

Lh,= a 

where 

L9 = [X m g - X aCcn + |I Zmg] CLa cn g (a w Eup + 'cn ' 'w) 
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The canard angle of attack is shown in Figure 4, where the subscript "up" indicates the 
effects of upwash due to the wing. 

The case for which there exists only a canard on the aircraft, a slightly different 
equation is derived due to the different directions the moments act about the point of 
ground contact. The C.G. of the aircraft is now put forward of the wing AC, which gives 
a negative moment created by the wing lift about the point of contact. Since the canard is 
placed forward of the C.G., the moment created by its lift needs to have a positive value to 
balance the moment equation. Equation 1 1 is derived from the summation of moments 



acn = aw + £up + icn - iw 
Figure 4: Forward Canard Angle of Attack 


about the point of ground contact with a forward canard and any thrust offset. 

M mg = - C Mwbti Sc ' [ Z mg" Z T 1 ] T l ' [ Z mg ' Z T^ T 2 * [ X mg X cg] W 

W - r It 

- [ X ac^- X mg] L wb + [ Z mg- Z t^ D + J V Z ™g + X mg " X ac cr | L cn 
Combining Equation 1 1 with Equations 4, 5, and 6, gives a solution for the canard lift 

- Lj - L 2 + L 4 + L 6 + L 7 + L 8 


(12) 


where Lj through L 6 , and L g are defined in the same terms as used in Equation 7, and 
where Ly and A are defined as 
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L7 — 


Xac wb ~ X mg - ^ Z mg , 


•'wb 


A = 


^ mg " ^ac cn + M- ^mg 


Finally, using the classic equation for the canard lift, the canard surface area is 
calculated using 


Sen = 


(13) 


Cloi g (a w + e up + icn ■ *w + 1 Sj q 


The above set of equations are used to determine the size of the horizontal control 
surface used for longitudinal control. These equations work for either an aft mounted 
horizontal tail, or a forward mounted canard. The next procedure is developed to insure 
that the horizontal stabilizer and control surface are large enough to permit an acceptable 
amount of C.G. travel while maintaining aircraft stability. 


CHAPTER 4 

Forward and Aft Center of Gravity Limits 

Using the preceding methods, it is possible to determine the minimum horizontal 
control surface size for take off rotation, and the C.G. range for the specified mission. 
Longitudinal controllability is now determined by calculating the forward and alt C.G. 
limits. These limits are then compared to the C.G. range determined. Then if necessary, 
the control surface size is increased to encompass the C.G. range within the forward and 
aft limits. 

First the forward C.G. limit is determined. This is calculated tor the worst possible 
case of flight, that is with the aircraft in the landing configuration. This includes effects 
due to high-lift at low-speed while in ground effect. Using the minimum horizontal 
control surface area calculated earlier, the equation for the summation of moments about 
the aircraft's C.G. is determined from Figure 5. This results in Equation 14. The forces 
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Figure 5: Summation of Moments about the Center of Gravity 


acting on the aircraft in the landing configuration are: the maximum wing lift, wing-body 
pitching moment with flaps extended, aircraft drag, horizontal control surface lift, and 
forces due to thrust offset. The wing lift is positioned at the AC. The wing-bodv 
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pitching moment is placed about the C.G. As was done earlier, the aircraft's drag is 
placed along the centerline of the aircraft (this force then falls out of the equation when 
assuming that the C.G. lies along the centerline). The horizontal control surface lift is 
placed at its respective AC, and the thrust offset forces are positioned according to the 
designer's inputs. 
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To find a solution for the forward C.G. position, the moment coefficient about the 
C.G. is set equal to zero, and Equation 14 is rearranged for determination ot X cg . The 

final solution is shown in Equation 15. By setting the moment coefficient to zero, a 
solution is calculated for the most forward C.G. position where the horizontal control 
surface is just able to maintain the aircraft in a level attitude with maximum elevator 
deflection. 


Xegforc = 
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In Equation 15, the C L is the maximum lift coefficient, and 5 c is the maximum elevator 
deflection. 

Once the forward C.G. limit is determined, it is compared to the different C.G. 
positions calculated for each of the mission phases of the aircraft. If any of these positions 
fall forward of the calculated limit, the horizontal control surface size is increased an 
incremental amount (5 square feet), and the forward limit is recalculated. This process is 
repeated until all the mission C.G. positions fall aft of the forward C.G. limit. 

The second comparison for the horizontal control surface size is the calculation ot the 
aft C.G. limit. The aft C.G. limit is that point where the moments about the C.G. no 
longer change with angle of attack. This is written as 

^Cm cg _ q ( 16 ) 

da 

This aft limit is also known as the neutral point (N 0 ) of the aircraft. If the C.G. moves 
behind this point, the change in moment with respect to angle ot attack becomes negative 
which makes the aircraft unstable. The N 0 of the aircraft is determined starting with the 

equation for the moment about the aircraft's C.G., Equation 14. Differentiating this 
equation with respect to angle of attack results in Equation 17. This includes contributions 
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( 17 ) 
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of the wing, horizontal tail, and forward canard. The fuselage and thrust vector 
contributions are small with respect to the other effects, and are therefore removed from 
the differentiation. 

The neutral point is determined using the equation derived from the following 
procedure. Equating C^ a to zero, and rearranging Equation 17 in terms of X cg results in 
Equation 18. In this equation, the N 0 is given in percent wing MAC. In order to compare 
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it to the C.G. range determined earlier, the aft C.G. limit must be calculated in units of 
length from the nose of the aircraft. The equation used for this calculation is shown in 
Equation 19. 

Xcg all = X LE + MAC N o ( 


Once the aft C.G. limit is determined, it is compared to the C.G. range. If any the 
C.G. positions do not fall forward of the aft limit, then the horizontal control surface is 
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increased in size an incremental amount (5 sq. ft.). The C.G. limit is reevaluated, and 
compared to the C.G. positions. This is repeated until the forward and alt limits 
encompass all the C.G. positions, and the aircraft becomes stable throughout all phases ot 

the mission. 

The horizontal control surface size is now the minimum size necessary to maintain 
aircraft control during take off rotation and low speed landing flight. It is also large 
enough to insure controllability, since the module assures the C.G. will remain in limits. 
The final step in sizing the horizontal control surface is to determine it it is large enough to 
insure aircraft stability during flight. 



CHAPTER 5 


Aircraft Stability 

In order to determine if the conceptual aircraft design is statically stable in flight, the 
pitching moment curve slope (dC Mcg /dC L ) or the static margin (SM) must be determined. 

The aircraft is stable if it has a negative value for the dC M /dC L , the larger the negative 

value, the more stable the aircraft. The equation used in solving for the dC M /dC L is 
developed from a summation of moments about the C.G., and is nondimensionalized by 
dividing by the dynamic pressure, the reference wing area, and the reference chord length. 
(Note that the drag does not appear since it is placed at the C.G. of the aircraft due to the 
limitations of the ACSYNT program in its calculations of total aircraft drag.) 
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which is modified by taking the derivative with respect to the lift coefficient. 
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Each section in the equation represents the contribution of the different aircraft components 
such as the wing, fuselage, aft and forward horizontal control surface, and engine thrust. 
The contribution of each component is calculated separately, and the results summed to 
determine the total aircraft stability coefficient. 

The wing contribution is solved using 4 


dCm = Xc g - Xac 
dC L wing v 


( 22 ) 


This solution assumes small angles of attack. The contribution of the wing is stabilizing 
when the C.G. is forward of the AC, and destabilizing when it is aft. 

The contribution of the fuselage and nacelles, dC Mwb /dC L , is estimated using the 

equation 5 6 

dC M wh Kr BDMAX 2 BODL (23) 

dCL S c Cloi 


where the fuselage stability coefficient K f , is determined from Figure 6. K f is given as a 

function of the position of the wing root quarter-chord. 

The contribution of an aft mounted horizontal control surface is solved using'’ 

dCm _ -_ML Sht Xacm ' ' X £S. TltlT ( i , de| (24 

dCL ht a S c da 


4 Perkins, C.D. and Hage, R.E., Airplane Pe rformance. Stability and Contr ol, J. Wiley 

and Sons, 1949. 216-218 

5 ibid., 229. 

6 ibid., 219-220. 
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which includes effects due to downwash of the wing on the tail. In this solution, 
stick-free effects are included by multiplying the horizontal control surface contribution by 

a factor called the free elevator factor (F e ), determination ot which is discussed in the 

chapter on "Additional Aerodynamic Calculations". 

The effect of a forward canard on the stability of the aircraft is calculated using 

dCm = acN Scn *CK n rN [ 1 - \ ( 25 > 

dCL cn a S c dcx ' U p 

The upwash of the wing on the canard, and the lift curve slope of the canard are calculated 
in the section "Additional Aerodynamic Calculations". The stick-lree effects are included 
in the canard calculations only if there is no aft horizontal tail on the aircraft. It both a 
canard and an aft mounted tail exist, then it is assumed the canard is used for trim only. 
The stick-free effects for the canard are calculated in the same manner as the stick-free 
effects for the horizontal tail. 



vs. 

Position of Wing Root Quarter-Chord (C iW4) 
(Reproduced from Reference 4) 
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The contribution of engine thrust is calculated using Equation 26, which assumes that 
the lift equals the weight of the aircraft. 

dC M -n, = ( TiZti+T 2 Zt2 \ n (26 

dC L ' W c ' 


Where the contribution of the thrust is due to offset of the thrust vectors from the aircraft 
centerline. This solution allows the designer to include effects for multi-engine aircraft 
where not all engines are located the same height from the centerline, and engines with 
more then one exhaust nozzle creating multiple thrust vectors. 

A second measure of aircraft stability is the SM. This is simply the distance the N () is 
aft of the C.G. at any given time during flight. It is an indication of how easy or hard it is 

for the aircraft to rotate about the lateral axis. The larger the distance between the N () and 
the C.G., the larger the SM, and the harder it is for the aircraft to be rotated. Conversely, 
the smaller the SM, the easier it is to rotate the aircraft. The SM is calculated from 


S M = 


N 0 - 


Xcg - XlE 

c 


(27) 


which gives the distance of the C.G. in front of the N 0 in percent MAC. Accepted values 

for the SM for stable aircraft range from 10% MAC for a transport, to 5 % MAC for a 
fighter (from reference 5). All the aircraft stability and control parameters have now been 
determined. There are however a few aerodynamic effects that influence the aircraft which 


need to be determined. 


CHAPTER 6 


Aerodynamic Effects 

In calculating the size of the horizontal control surface needed to maintain stability, 
certain aerodynamic characteristics that affect the solution need to be determined. Three 
important characteristics are included in the stability module. These include the downwash 
of the wing on an aft horizontal tail, the upwash of the wing on a forward canard, and 
ground effects. 

Following a method outlined in reference 1, the downwash of the wing on an aft 
mounted horizontal tail is calculated 7 . This method gives the downwash angle at the tail as 

a function of the effective wing aspect ratio (AR cff ), the effective wing span (b cff ), and the 
tail height above or below the trailing wing vortex. This method assumes a large wing 
span to horizontal tail span ratio (b/b h[ > 1.5), wing trailing edge vortex separation, and 
subsonic flow. 

First AR cff and b, ff are determined. These are determined from Figure 7, dependent 

on the angle-of-attack and geometry of the wing. Once the AR cff and the b cff have been 
determined, the downwash angle is calculated. Equations 28, 29, and 30 are used to solve 
for the span of the vortex core (b v ) at the horizontal tail quarter-chord. 


b v = b 
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2 Hoak, D.E. et al; USAF Stability and Control Datcom ; Wright Patterson AFB Ohio, 
45433; Revised 1970; sect. 4.4.1 
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where L cff is the distance from the wing tip trailing edge to the horizontal tail quarter-chord 
and is a dummy variable. The following are used in solving for b v . 


and 



0.78 + 0.10 (A. - 0.4) + 0.003 A c/4 


b c ff 


(29) 
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0.56 AR cf f 

cZ 


(30) 


Equation 31 is used to solve for the height of the horizontal tail above or below the trailing 
wing vortex core (a). This is a function of the tail height above the wing chord line (h HT ), 
and the effective distance from the wing vortex separation to horizontal tail. 


a - h^u - L cff 

Finally, the ratio of the downwash at 
calculated using 


1 rc AR cff 


the tail to the downwash at 


(31) 

the vortex core is 



(32) 


where the downwash at the vortex core is determined from 


1 .62 C L 
7t AR 


(33) 


where the AR for this equation is the actual wing AR, and not effective. 

The upwash effects of the wing acting on a forward canard are determined as a 

function of wing AR, wing root chord (C rool ), and the distance the canard is ahead of the 
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wing quarter-chord. The upwash effects are determined using Figure 8 8 . This solution is 
accurate only for Mach numbers less then one, and as such currently limits the analysis to 
subsonic missions. 



Figure 8: Upwash Effects in Front of the Wing(l - d e/da) 

vs. 

Forward Posiuon in Percent Wing Quarter-Chord (x/c root) 
as a Function of Wing Aspect Ratio (AR) 


Since principal sizing of the horizontal control surface is computed at takeott or 
landing, it is important to include ground effects. Ground effects can adversely at feet the 
horizontal control surface size by decreasing the downwash or upwash angles on the 
horizontal control surface. This decreases the overall angle of attack of the surface, which 
then increases the control surface size needed to generate the required amount ot lilt. 
There are three areas which are influenced by ground effect. The first two are the lift ot 


8 McCormick, B.W., Aerodynamics. Aeronautics, and Flight Dyn amics, J. Wiley and 
Sons, 1979. pg. 520. 
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the wing and the lift of the horizontal control surface. The third is the downwash angle of 
the wing on the horizontal control surface. 

As discussed earlier, when sizing the horizontal control surface, it is necessary to 
calculate the lift produced by the wing and the lift required of the horizontal control 
surface. Once these have been determined, it is a simple matter of modifying them to 
include ground effects. Ground effects increase the closer the aircraft is to the ground. To 
get maximum effect of the ground on the horizontal control surface size, the height of the 
aircraft above the ground is reduced to a minimum. The minimum possible height occurs 
when the aircraft is resting on the ground with landing gear extended. This height is used 
along with Figure 9 to calculate the ratio of lift-curve slopes in and out of ground effect 9 . 
Once the change in the lift curve slope is determined, the wing lift in ground effect is 
calculated using 

L w b=C Lag «w4S (34) 


the horizontal tail lift is determined using 
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(35) 


and the canard lift is found using 

L cn = C Lacn e up + *cn - *S c )q S cn (36) 

In determining how ground effects influence the downwash angle, the method used 
calculates the change in downwash as a function of several wing and horizontal tail 


9 Perkins, C.D. and Hage, R.E.; Airplane Perf ormance. Stability and Control; 
J. Wiley and Sons, 1949; 257. 
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geometric and aerodynamic parameters 10 . These include effective wing span, wing aspect 
ratio, wing taper ratio, wing height above ground, and horizontal tail height above 
ground. These also include the wing lift coefficient, and the change in wing lift coefficient 
due to flaps. The effective span of the wing is determined as a function of the wing lilt, 



_d£_ 

b/2 


Figure 9: Ratio of Lift Curve Slopes in and 
out of Ground Effects (a^a) 
vs. 

Height above Ground in Semi-Spans (dg/b/2) 
as a Function of Aspect Ratio (AR) 


the change in wing lift due to flaps, and the effective wing and flap span ratios, as 
calculated using Equation 37. Once the effective span is calculated, the change in 
downwash is determined using Equation 38, where H is the wing height above the 


10 Hoak, D.E. et al; 1ISAF Stability and Control Datcom ; Wright Patterson AFB Ohio, 
45433; Revised 1970; sect. 4.7.1 
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ground, and H m is the height of the horizontal tail above ground, 
change in downwash angle is good for Mach numbers less than one. 
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The ratios b w '/b, and b r 7 b w ' are determined using Figures 10 and 
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CHAPTER 7 
Module Verification 


The stability module was validated by comparing computer generated data of two 
types of aircraft to data available on the production aircraft. Computer models of a 
General Dynamics F-16A and a Boeing 727-200 where generated for comparison using 
the ACSYNT program. These computer models matched various parameters of the 
production aircraft which included geometry, aerodynamics, and weights. A comparison 
was then made between the horizontal control surface size for the production aircraft 

versus the control surface size calculated using the stability module. 

The method used to solve for the horizontal control surface size consists of three main 
steps. First, the different aircraft components are placed along the fuselage. It is 
important to place the major components of the aircraft in their specific locations. These 
include the wing, horizontal control surface, engine(s), fuel in the fuselage, fuel in the 
wing, and the main landing gear. Once the components are in place, the C.G. position 
and range for the actual aircraft is matched in the computer model, allowing tor a shift in 
C.G. from the fully forward to fully aft positions. This matching necessitates the moving 
of less important components to different positions, or recalculating the weights of some 
components. After determining the C.G., the ACs of the wing and horizontal control 
surface are compared to assure the lift forces are acting at the correct positions. 

Once the C.G. position, C.G. range, and the ACs correspond to their respective 
positions on the actual aircraft, the horizontal control surface sizes are compared. The 
results, which are presented in Tables 1 , and 2, show excellent agreement between the 
actual and computer aircraft. Since more data was available for the F-16A, the solution ot 
the computer module matched best with the production aircraft. 
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To validate the module for an aircraft with a forward mounted canard, the geometry 
and weight of the F- 16A model was modified to match the dimensions of a Saab Viggen. 
The comparison for this model is shown in Table 3. This showed good accuracy, with an 

error in the canard size of only 3.91 percent. 

An important factor in these tables is the position of the main landing gear. The 
stability module uses rotation about the main landing gear to determine the size of the 
horizontal control surface. It was found that the position of the main landing gear at take 
off had a large effect on the calculated size. Using the F-16A as an example. Figure 12 
shows that small changes in main landing gear position can effect the calculated control 
surface size. Figure 13 shows how changes in main landing gear position atfect the 
forward and aft C.G. limits. Note in both of these figures the tail size and the C.G. limits 
stabilize when the main gear fall at or forward of the 60 percent fuselage length. This is a 
result of the stability module forcing the main gear aft to satisfy the 15 degree angle 
requirement between the mnin gear and the C.G. 
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Table 1 : F-16A Production and Computer Model Comparison 



Production 

Computer Model 

% Difference 

Horizontal 
Tail Size 

49.0 

48.18 

1.67 

Main Gear 
Position 

29.0 

28.69 

1.07 

MAC 

11.3 

11.20 

0.88 

Quarter-chord 

Position 

25.6 

25.55 

0.20 

Half-chord 

Position 

28.4 

28.35 

0.18 

Forward C.G. 

* 

26.10 


Aft C.G. 

* 

27.21 


Quarter-chord of 
Horizontal Tail 

41.5 

41.32 

0.43 


Note: All distances are in feet from the nose of the aircraft, except the Mean 
Aerodynamic Chord is in feet from the leading to trailing edges. 

Horizontal tail size is in feet squared. 


*: The author was unable to obtain information on the 
forward or aft C.G. limits of the Production aircraft 
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Table 2: Boeing 727-200 Production and Computer Model Comparison 



Production 

Computer Model 

% Difference 

Horizontal 
Tail Size 

376.0 

379.64 

0.96 

Main Gear 
Position 

66.18 

68.47 

3.34 

MAC 

15.5 

13.98 

9.81 

Quarter-chord 

Position 

63.3 

62.84 

0.73 

Half-chord 

Position 

67.2 

66.34 

1.28 

Forward C.G. 

* 

62.14 


Aft C.G. 

* 

65.19 


Quarter-chord of 
Horizontal Tail 

113.7 

113.44 

0.23 


Note: All distances are in feet from the nose of the aircraft, except the Mean 
Aerodynamic Chord is in feet from the leading to trailing edges. 

Horizontal tail size is in feet squared. 


*: The author was unable to obtain information on the 
forward or aft C.G. limits of the Production aircraft 
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Table 3: Saab JA-37 Viggen Production and Computer Model Comparison 



Production 

Computer Model 

% Difference 

Horizontal 

CanardSize 

66.7 

64.09 

3.91 

Main Gear 

34.8 

36.14 

3.71 

Position 



MAC 

17.15 

17.94 

4.40 

Quarter-chord 

35.3 

32.64 

7.54 

Position 


Half-chord 

39.64 

37.13 

6.33 

Position 



Forward C.G. 

* 

31.40 


Aft C.G. 

* 

36.05 


Quarter-chord of 
Horizontal Tail 

- 

- 

- 


Note: All distances are in feet from the nose of the aircraft, except the Mean 
Aerodynamic Chord is in feet from the leading to trailing edges. 

Horizontal cananrd size is in feet squared. 


*: The author was unable to obtain information on the 
forward or aft C.G. limits of the Production aircraft 




Figure 12: Horizontal Tail Size (Sht) 
vs. 

Main Landing Gear Position (Xmg) 




CHAPTER 8 


Vectored Thrust Analysis 

The stability and control module includes a subroutine which allows for the vectoring 
of thrust during landing. This gives the designer the ability to evaluate aircraft with 
hovering and vertical landing capabilities, using the stability equations to create stable 
designs. 

This subroutine is used to calculate one of three different parameters tor a stable 
aircraft in high-lift, low-speed transitioning flight. The three parameters include the 
forward thrust vector position (X T1 ), the aft thrust vector angle (Y 2 ). and the thrust 

split between the forward thrust vector and the total thrust available (TSPLI T). The 
stability module is used to solve for these three parameters because they have the greatest 
impact on the aircraft's design. The forward thrust vector position has a great effect on the 
internal arrangement of the aircraft. It represents a large lift producing system, either 
ducting from the main engine, or a separate auxiliary lift engine. It therefore requires a 
large amount of internal volume, limiting the placement of other aircraft components. The 
aft thrust vector angle and the thrust split determine the amount of bleed air that is removed 
from the main engine. This has a direct effect on the size and thrust of the engine. The 
more bleed air removed from the engine, the larger it needs to be. 

The forward thrust angle is not considered important since the most effective thrust 
angle is at 90 degrees to the ground. The aft thrust vector position is not considered 
important in the stability module because it is more economical both in design and weight 
of the engine to put the aft thrust at the rear of the aircraft (as in conventional aircraft). 

The general layout for the thrust vector angles and distances is shown in Figure 14, 
which shows distances for both the C.G. and the main gear. During take off and cruise 
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flight, it is assumed the thrust vectors point in the aft direction, as with a conventional 
aircraft. This allows the designer to include any thrust offset due to the positions of the 
engines. During landing, the thrust nozzles are "rotated" and the thrust forces are applied 
in both the horizontal and vertical directions. 



Figure 14: Thrust Vector Components 
and Positions 

Summation of the thrust and aerodynamic forces about the C.G. result in the familiar 
equation for C Mcg (Equation 14). This now includes the thrust forces shown in Figure 
14, the final result shown in Equation 39. 
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In order to maintain stability, the summation of moments about the C.G. is set equal to 
zero. Solutions can now be found for the forward thrust vector position, the aft thrust 

vector angle, and the thrust split. 

Solving for X T1 is a straight forward algebraic rearrangement of the moment equation 

given above. Equation 40 is used to solve for the forward thrust vector position in units 
of length from the nose of the aircraft. This solution is for a stable aircraft in low speed 
flight transitioning from forward flight to hover. 

_ SUMVT q S c x (4 0) 

Aji - ' eg 


where SUMVT is given by 
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By definition, the thrust split (TSPLIT) is the ratio of the forward thrust over the total 
thrust. A solution for the thrust split can be determined by using relationships between the 
forward thrust, the aft thrust, the total thrust, and the thrust split. The first relationship 
comes from the definition of the thrust split, as shown in Equation 41. The second 


TSPLIT = 


Ti 

THRUST 


(41) 


relationship comes from the assumption that the aft thrust equals the total thrust minus the 
forward thrust (this ignores any frictional and heat losses in the forward thrust). 
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T 2 = THRUST - Tj = THRUST (1 - TSPLIT) 

Equations 41 and 42 are substituted into Equation 39, which is rearranged to solve for 
TSPLIT. The final solution is shown in Equation 43. 


TSPLIT = 


SUMVT Y sin ?2 7 cos "h 
THRUST ' XT2 q Sc T2 q S c 

A 


(43) 


where 


SUMVT = C L 


(Xcg-Xaj „ 

- c M wh - L L« a Hr s 


S ffp (Xac 1(r X cg ) 


HT 


+ C, a 

L aCN 


S CN ( 


x ac -x ce 

ac CM c & j 


CN S 


^CN 


and 


cosy 2 siny 2 siny, cosy, 

A = ZT2 "qS^ " Xt2 qTc ' Xti q S c " T1 q S c 

Solving for the aft thrust vector angle (y 2 ) is more involved because it appears in the 

moment equation twice, once in a sine function, and once in a cosine function. Since a 
direct solution can not be found, an iterative approach is used. The boundaries of the 

solution are known to be 0 degrees and 90 degrees, and are therefore as starting points. 

The moment equation, Equation 39, is solved for C Mcg with y 2 set to 0 degrees, and with 
y 2 set to 90 degrees. The two solutions are compared, and the one that has the largest 

magnitude for C Mcg is reduced to halfway between the two boundaries. The moment 
equation is recalculated for the remaining boundary and for the new boundary value. 
Again, the solution with the largest magnitude for C Mcg is reduced to halfway between the 
two boundaries. This series of computations and comparisons continues until the 
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magnitude for the moment equation approaches zero. The corresponding Y 2 is a solution 

for a stable aircraft design for transition from horizontal flight to hover. 

The vectored thrust subroutine was used to compare the forward thrust vector 
positions to the thrust split, and the aft thrust vector angle. The results of this can be seen 

in Figure 15. 



Xti (feet) 



Figure 15: Forward Thrust Vector Position (Xti) 
vs. 

Thrust Split (TSPLIT) for Various 
Aft Thrust Angles ( 72 ) 


* This graph is for an F-16A modified with 

with a forward nozzle set at y= 90 degrees and 
with the aft nozzle set at a position of 
0.95 percent of body length. 




CHAPTER 9 

Additional Aerodynamic Calculations 


In the ACSYNT design program, most of the aerodynamic parameters needed by the 
stability module are solved for in the aerodynamics module. However, some variables are 
not solved for and others are currently impossible to transfer between the two different 
modules. These additional variables are calculated inside the stability module using the 

methods discussed, as taken from the specified references. 

The aerodynamic center (AC) of the wing is the first of the aerodynamic variables to 
be determined. In the aerodynamics module of ACSYNT, the AC is assumed to be the 
quarter chord point of the MAC of the wing, dependent only on the geometry of the wing 
and not the Mach number of the aircraft. This is not an accurate representation when 
designing high speed fighter aircraft. The aerodynamics module does, however, 
determine the change in the lift-curve slope of the wing with a change in Mach number. 
This information can be used to determine the change in AC with the change in Mach 
number with a fair degree of accuracy. 

Figure 16 shows the forces and distances used in the calculation of the AC . The 
summation of moments about point A results in Equation 44. This equation is 
nondimensionalized by dividing by the dynamic pressure, the wing chord, and the wing 
area. This equation is then reduced by taking the derivitive with respect to the angle-of- 

b 

M a = q f c 2 C Mac dy - q 


1 1 McCormick, B.W., Aerodynami cs. Aeronautics, and Flight Dynami cs, J. Wiley and 
Sons, 1979. pg. 484 - 485. 


L 


c C, y tan A dy 


(44) 


44 


45 



Figure 16: 

Calculation of Wing Aerodynamic Center 
by Moments about Root Quarter-Chord, point A. 


attack. Noting that by definition the moment coefficient about the AC (CM ac ) does not 
change with angle-of-attack, it is therefore aero. This reduced form of this equation is 
shown in Equation 45. 


dC M A 

d a 



y tan A dy 


2 


(45) 


Defining X A to be the distance from the point A to the AC, the the moment about the 
AC (M ac ) is defined as the summation of the moment about point A (M A ), and the lift 
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force multiplied by the distance X A . This is seen in Equation 46. 

M ac = M a + LX a 


( 46 ) 


This equation can be nondimensionalized by dividing by the dynamic pressure, the wing 
chord, and the wing area. This can also be reduced by taking the derivative with respect to 
angle-of-attack. This results in Equation 47, again noting that the moment about the AC 
does not change with change in angle-of-attack. 

d Cm * _ ^ (47) 

— '-La r 


The distance X A is determined using the equation that results from the combination of 
Equations 45 and 47, as shown in Equation 48. 


X A = 


C^S 




c Cj a y tan A dy 


(48) 


This general equation is simplified by assuming a constant lift curve slope along the 
span of the wing (thereby removing it from the integral), and by assuming a linearly 
tapered wing (thereby removing the integral completely). This simplified form is shown 
in Equation 49. 


X 


ac 


1+2*. 1 /b . 

tttM* 


'la 


'La 



^root 

“4~j 


(49) 


This equation is used to calculate the position of the AC from the point A (the 
quarter-chord point of the wing center), and accounts for wing sweep and taper ratio. The 
final step is to account for the effects due to Mach number. 
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The only aerodynamic characteristic in this equation is the ratio of the 
two-dimensional to three-dimensional lift curve slopes. It is therefore necessary to 
determine the change in the lift curve slope with respect to Mach number. This works 
well with the ACSYNT design program since the only variable that ACSYNT determines 
Mach effects on is the three-dimensional lift curve slope of the wing. The stability module 
therefore needs only to be used to determine the ratio of the two-dimensional to 
three-dimensional lift curve slopes to include Mach effects on the AC. 

The ratio of two-dimensional to three-dimensional lift curve slopes is determined by 

rearranging Equation 48 . 

r = El? (50) 

^ 57.29 C la / 

1 /e 7t AR 

The AC is now determined for any geometric shape, and at varying Mach numbers. 

It was found that this method is good for Mach numbers up to Mach 1.2, at which 
point the AC changes at a rapid rate. Comparisons between the method used in the 
stability module and a graphical method discussed in Reference 1, are shown on Figure 
17. This shows the AC shift for a Boeing 727-200 in subsonic flight, and a F-16A in 
subsonic and supersonic flight. The equations used to determine the AC tend to give a 

higher value, averaging 4 % for the 727-200. 

The second aerodynamic parameter not determined in the aerodynamics module of 
ACSYNT is the three-dimensional lift curve slope of the canard. A solution is found first 
using the Helmbold equation 13 , Equation 51. This gives the three-dimensional lift curve 


12 Perkins, C.D. and Hage, R.E., Airplane Pe rformance. Stability and CbntmL 

J. Wiley and Sons, 1949. pg. 220. 

13 McCormick, B.W., Aerodynamics. Aeronautics, and Flight Dyna mics, J. Wiley and 

Sons, 1979. pg. 137. 
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— Numerical Results 
Graphical Results 



B727-200 


— Numerical Results 
Graphical Results 



F-16A 


Figure 17: Numerical and Graphical Aerodynamic 
& Center Comparisons 

(Graph Reproduced From Reference 1) 
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slope as a function of the two-dimensional lift curve slope and the AR of the canard. 
Equation 51 is modified to include the effects of canard sweep angle and Mach number 
through modification of the two-dimensional lift curve slope 14 . Sweep angle effects are 
included by the multiplication of the two-dimensional lift curve slope and the cosine of the 
sweep angle. The Mach number effects are included through division of the 
two-dimensional lift curve slope by the Pradtl-Glauert compressibility factor. The 
compressibility factor depends on the subsonic or supersonic flow. The linal 
three-dimensional lift curve slope equation is 


Clu - 


C la AR 



- (AR M^ 2 


(52) 


This equation is simplified somewhat in the stability module by assuming the 
two-dimensional lift curve slope of the canard equals two-pi. 

Also in the stability module, the contribution of the horizontal control surface is 
modified to include the effects of a stick-free condition. This is determined by calculating 
the hinge moment parameters of the elevator. The method is used to modify the 

contribution of the horizontal control surface to the dC M /dC L , and the equations for the 


control surface lift. This is done by multiplying with the free elevator factor (F c ) 15 , 
which is calculated from 


i d8 bj 

= 1 - x = 1 - Xt— 

da b 2 


(53) 


14 


ibid., 283 - 284. 
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Where b, and b 2 are determined from the multiplication of the coefficients found in 
Figures 18, 19, 20, and 21. These are given as functions of the elevator span, t/c, AR, 
and balance ratio (BR), each specified of the designer. 

bj = -0.55 k l(cc/c) k l(t/c) k l(BR) k l(1/A) (54) 

b 2 = -0.89 k 2(Cc/c) k 2(t/c) k 2(BR) k 2(1/A) (55) 

The control effectiveness factor (X), is determined from Figure 22 as a function of elevator 
to stabilizer chord ratio. 


15 ibid., 495 - 508. 
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Figure 22: Control Surface Effeciency Factor (t) 
vs. 


Percent Elevator Chord (ce/c) 



CHAPTER 10 

Conclusions and Recommendations 


The use of a design program can enhance many areas in the conceptual aircraft design 
stage. NASA Ames Research Center has created a computer program, ACSYNT, which 
does calculations of aircraft geometry, aerodynamics, propulsion, mission performance, 
and weights. This allows designers to examine a wide range of designs in a relatively 

short period of time, to create the best possible aircraft. 

A stability and control module was created to enhance the conceptual design program, 
ACSYNT. This module calculates the size of the horizontal control surface, the center of 
gravity for each mission phase, the forward and aft center of gravity limits, and the 
longitudinal stability of the design. The stability module calculates the control surface size 
needed for take off rotation, it examines the stability of the aircraft during the mission, and 
it determines if the tail size is large enough to maintain controllability during landing. If 
the tail size is not large enough to meet any of these goals, it is increased until stability and 

controllability are established. 

Comparisons between production aircraft and computer models show that the stability 
module accurately sizes the tail size for a range of aircraft types. Three aircraft, a 
B727-200 transport, and a F-16A fighter, and a JA-37 Viggen fighter with a forward 
mounted canard, were used to evaluate the stability module. In each case, the tail size 
determined by the module was an accurate representation of the actual aircraft. 

In addition to analyzing the horizontal control surface size and the stability of the 
aircraft, an analysis was done to evaluate vectored thrust applications on the design of a 
conceptual aircraft. This allows the designer to evaluate one of three different parameters 

associated with vectored thrust systems, while maintaining stability of the aircraft during 
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landing. These include the forward thrust vector position, the aft thrust angle, and the 
thrust split. These parameters were determined to be the most important in the design of 
the aircraft. The forward thrust position affecting the internal layout of the aircraft, and the 
aft thrust vector angle and the thrust split affecting the required thrust of the engine. 

Recommendations for improvement on this module include the following. 

- Calculations of aerodynamic center at higher Mach numbers than what is currently 

being used. 

- Include calculations to determine the downwash and upwash effects of the wing at 
supersonic Mach numbers. 

- Calculations of ground effects on the upwash of the wing. 

- Modification of wing lift to include interaction between the wing and the canard 
vortices. 

- In the Vectored thrust analysis, include the ability to augment the forward thrust 
with an afterburner, duct burner, or ejector. 

- In the Vectored thrust analysis, include changes in lift and pitching moment due to 

jet interaction. 
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APPENDIX A 


Module Inputs 


The following is a list and description of inputs used by the stability and control 
module. All the input variables have default values, so not all of the variables need to be 
input. The user inputs these in the namelist 'STABIN' and select.ng the module number 5 

in the COPES inputs. 

Format for the 'STABIN' namelist includes a title line, a maximum of 80 characters 
long, followed on the next line by 'SSTABIN' and then the desired inputs. Finally, at the 
end of the inputs, a 'SEND' statement is needed to tell ACSYNT the input has been 

completed. 

Example: 

******* Stability and Control Inputs, F-16A Falcon ******* 

SSTABIN ETACN=.90, XT1=0.40, 

XFCREW=0.25, XFINST=0.21, 


SEND 


Name 

Default 

Description 

— Keai iormai 


AWTOT 

2.00 

Angle of attack of the wing at take off rotation 

CBALHT 

0.00 

Percent chord of the longitudinal control 



surface ahead of the hinge line. 

CELV 

0.30 

Percent chord of the horizontal stabilizer 



that is the movable surface. 

ETACN 

0.90 

Canard efficiency. 

ETAHT 

0.90 

Horizontal tail efficiency. 


(degrees) 
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ELVDEF 

-10.00 

Amount of elevator deflection needed for 
takeoff rotation, (degrees) 

ELVDMAX 

-15.00 

Maximum amount of elevator deflection needed for 
controllability at low speed, (degrees) 

GAMAT1 

0.00 

Forward thrust nozzle rotation angle from 
the horizontal, (degrees) 

GAMAT2 

0.00 

Aft thrust nozzle rotation angle from the 
horizontal, (degrees) 

ICN 

0.00 

Canard incidence angle, (degrees) 

IHT 

0.00 

Horizontal tail incidence angle, (degrees) 

SCN1 

0.0 

Area of canard when both a canard and an aft tail 
are being used (feet 2 ). 

SPANCN1 

0.0 

Span of the canard when both a canard and an alt 
tail are being used (feet). 

TRCN1 

0.0 

Taper ratio of the canard when both a canard and 
an aft tail are being used. 

TSPLIT 

0.00 

Percent thrust split, front- to-rear. 

TZCHT 

0.10 

Elevator thickness-to-chord ratio. 

VROT 

0.30 

Mach number at take off rotation. 

WLGFRT 

0.30 

Weight ratio of front-to-rear landing gear 
for C.G. calculations. 

XFAMMUN 0.30 

C.G. position of the ammunition in percent BODL (fuselage 
length). 

XFAPU 

0.95 

C.G. position of A.P.U. in percent BODL. 

XFBB2 

0.50 

C.G. position of advanced weapons in percent BODL. 

XFBOD 

0.50 

C.G. position of body in percent BODL. 
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XFBOMB 

0.50 

C.G. position of bombs in percent BODL. 

XFCAN 

0.20 

C.G. position of the canard in percent of the Mean 
Aerodynamic Chord from the leading edge. 

XFCARGO 

0.50 

C.G. position of cargo in percent BODL. 

XFCREW 

0.10 

C.G. position of crew in percent BODL. 

XFELT 

0.05 

C.G. position of Avionics in percent BODL. 

XFENG 

0.50 

C.G. position of the engine in percent of PODL (engine 
pod length). 

XFEP 

0.50 

C.G. position of electrical system in percent BODL. 

XFEX 

0.50 

C.G. position of fuselage external tank in percent BODL. 

XFFFUS 

0.50 

C.G. position of fuselage fuel in percent BODL. 

XFFS 

0.60 

C.G. position of fuel system in percent BODL. 

XFFUR 

0.50 

C.G. position of furnishings in percent BODL. 

XFFWG 

0.50 

C.G. position of the fuel in the wing in percent BODL. 

XFHDP 

0.50 

C.G. position of the hydraulics and pneumatics 
in percent BODL. 

XFHT 

0.20 

C.G. position of the horizontal tail in percent of the Mean 
Aerodynamic Chord from the leading edge. 

XFINST 

0.15 

C.G. position of the instruments in percent BODL. 

XLECN 

0.0 

Leading edge position of canard when both a 
canard and an aft tail are used (feet). 

XFLIFTF 

0.50 

C.G. position of the lift-fan in percent BODL. 

XFLGFRT 

0.10 

C.G. position of the nose gear in percent BODL. 

XFLGR 

0.55 

C.G. position of the main gear in percent BODL. 

XFMISS 

0.50 

C.G. position of the missiles in percent BODL. 

XFNA 

0.98 

C.G. position of the nacelles in percent PODL. 
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XFPA 

0.55 

C.G. position of passenger accommodations 
in percent BODL. 

XFPAYL 

0.50 

C.G. position of any payload in percent BODL. 

XFPIV 

0.65 

C.G. position of the control surface pivots 
in percent BODL. 

XFSC 

0.70 

C.G. position of control surfaces in percent BODL. 

XFVT 

0.20 

C.G. position of the vertical tail in percent of the Mean 
Aerodynamic Chord from the leading edge. 

XFWG 

0.15 

C.G. position of the wing in percent of the Mean 
Aerodynamic Chord from the leading edge. 

XT1 

0.30 

Horizontal position of the forward thrust vector in percent 
BODL. 

XT2 

0.60 

Horizontal position of the aft thrust vector in percent BODL 

ZFD 

0.00 

Vertical distance of the aircraft centerline from the C.G. in 
percent BDMAX (maximum fuselage length). 

ZFMG 

2.00 

Vertical distance of the main gear from the C.G. in percent 
BDMAX. 

ZFT1 

0.00 

Vertical distance of the forward thrust vector from the C.G. 
in percent BDMAX. 

ZFT2 

0.00 

Vertical distance of the aft thrust vector from the C.G. in 
percent BDMAX. 

ZMU 

0.02 

Friction coefficient of the runway on the 
main gear. 

ZRTCN1 

0.0 

Height of the canard quarter-chord above the aircraft 
centerline (feet). 
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— Integer Format — 


IAS 

0 

Trajectory mission phase in which ammunition is used. 

IBS 

0 

Trajectory mission phase in which bombs are dropped. 

ICGPRT 

0 

Print flag to output aircraft component C.G. 
positions for C.G. analysis. 

IDBPRT 

0 

Print flag to output a greater amount of information during 
each mission phase, used for debugging purposes. 

IETANK 

0 

Used to specify various positions of any external tanks on 
the aircraft. 


0 

a single tank at the fuselage centerline. 


1 

Two tanks mounted on the sides of the fuselage. 


2 

Two tanks mounted on the wing quarter-chord. 


3 

Two wing-tip tanks. 


4 

One fuselage centerline tank and two tanks mounted at the 
wing quarter-chord. 


5 

One fuselage centerline tank and two wingtip tanks. 


6 

Two wing quarter-chord and two wingtip tanks. 

IMS 

0 

Trajectory phase in which missiles are used. 

IVECT 

0 

Used to specify which vectored thrust parameter is to be 
solved for. 


0 

No vectored thrust variable is to be solved for. 


1 

Calculate longitudinal position of the forward thrust vector. 


2 

Solve for the Forward over aft thrust ratio (Thrust Split). 


3 

Solve for the angle of the aft thrust vector from the 
horizontal axis. 

NPHASE 

1 

Total number of trajectory phases. 


APPENDIX B 
Sample Output 
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AIRCRAFT COMPONENT WEIGHTS AND POSITIONS 


COMPONENT 

WIEGHT 

POSITION 

AIRFRAME 

17022.544922 

42.939999 

WING 

17129.294922 

64.515587 

CANARD 

0.000000 

0.000000 

HT 

1581.655640 

114.010048 

VT 

1310.330078 

107.414925 

NOSE GEAR 

1223.618774 

11.300000 

MAIN GEAR 

2855.110352 

70.512001 

NACELLS 

1503.169800 

45.200001 

PIVOTS 

0.001953 

45.200001 

AIR COND. 

309.482300 

44.365093 

APU 

605.342590 

70.059998 

AVIONICS 

1441.599487 

1.130000 

ELECTRICAL 

1902.595215 

22.600000 

INSTURMENTS 

796.524597 

5.650000 

HYDRAULICS 

563.688660 

56.500000 

CONTROL SURF 

2981.909424 

56.500000 

PA 

0.000000 

45.200001 

FURNISHINGS 

450.000000 

36.160000 

ENGINES 

6954.871582 

91.914490 

FUEL SYS. 

1682.630249 

33.900002 

LIFT FAN 

0.000610 

56.500000 

PAYLOAD 

60.000000 

22.600000 

CREW 

510.000000 

11.300000 

CARGO 

0.000000 

33.900002 

AMMUNITION 

0.000000 

33.900002 

BB2 

0.000000 

56.500000 

BOMBS 

0.000000 

56.500000 

MISSILES 

0.000000 

53.110001 

FUS. FUEL 

0.000 

56.500000 

WING FUEL 

16347.398 

64.235992 


MISSION PHASE NUMBER 1 
XCG = 64.0515594 

XAC = 64.0994263 

DCMDCL = -0.0539999 

NEUTPT = 0.4538571 

STATIC MARGIN = 1 1 .7050486 

AFT CG = 65.6878738 

SUMMO = 5729981.000000 
DCMDCLT = -0.05395812 

DCMDCLF = 0.00338223 

DCMDCLJ = 0.00000000 

LWB = 33122.109375 LHT = 

SHT = 379.950439 SCN = 

XCGFRWD[ 1] = 61.9601 


SUMWT = 89458.882813 

DCMDCLW = -0.00342405 

DCMDCLCN = 0.00000000 


-9019.528320 

0.000000 


LCN = 0.0000000 


B727-200 Center of Gravity Output and Detailed Output 



STABILITY OUTPUT DATA 


TAIL SIZE = 

XCGFRWD = 

AFT CG = 

DATA AT END OF TAKE-OFF 
XAC = 

NEUTPT = 

XCG = 

CGBAR = 

DCMDCL = 

STATIC MARGIN = 
MISSION PHASE NUMBER 
XAC = 

NEUTPT = 

XCG = 

CGBAR = 

DCMDCL = 

STATIC MARGIN = 
MISSION PHASE NUMBER 
XAC = 

NEUTPT = 

XCG = 

CGBAR = 

DCMDCL = 

STATIC MARGIN = 
MISSION PHASE NUMBER 
XAC = 

NEUTPT = 

XCG = 

CGBAR = 

DCMDCL = 

STATIC MARGIN = 
MISSION PHASE NUMBER 
XAC = 

NEUTPT = 

XCG = 

CGBAR = 

DCMDCL = 

STATIC MARGIN = 
MISSION PHASE NUMBER 
XAC = 

NEUTPT = 


XCG = 

CGBAR = 
DCMDCL = 

STATIC MARGIN = 


379.9504395 sq. ft. 
62.3969193 ft. 
65.0767975 ft. 

— F 

63.1220779 ft. 
0.4101450 
64.0528641 ft. 
33.6899948 % chord 
0.0047093 
7.3245034 % chord 
1 

64.0994263 ft. 
0.4538571 
64.0515594 ft. 
33.6806602 % chord 
-0.0539999 
7.3338361 % chord 
2 

64.4846954 ft. 
0.4740317 
64.0514832 ft. 
33.6801147 % chord 
-0.0784065 
7.3343816% chord 

3 

64.1274567 ft. 

0.4552816 
64.0369263 ft. 
33.5759850 % chord 
-0.0568189 
7.4385114 % chord 

4 

64.5148773 ft. 
0.4756619 
64.0367661 ft. 
33.5748367 % chord 
-0.0814076 
7.4396572 % chord 

5 

64.2071304 ft. 
0.4593702 
64.0241013 ft. 
33.4842453 % chord 
-0.0627594 
7.5302525 % chord 


END STABILITY OUTPUT 


B727-200 Final Stability Output 



APPENDIX C 
Production Aircraft Layouts 
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Figure 26: B727-200 Aerodynamic Data 
(Reproduced from Reference 7) 
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Program Listing 
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SUBROUTINE STBLTY 


THIS SUBROUTINE IS USED IN CONJUCTION WITH NASA AMES AIRCRAFT DESIGN 
PROGRAM ’ACSYNT ’ . THIS PROGRAM CALCULATES THE SIZE OF THE HORIZONTAL 
CONTROL SURFACE FOR EITHER AN AFT TAIL OR A FORWARD CANARD^ IT CAN 
ALSO HANDLE AN AIRCRAFT WITH BOTH CONTROLS, BUT IT ONLY SIZES THE AFT 
TAIL THE SIZE OF THE CANARD IS AN INPUT FOR THE USER 
THIS SUBROUTINE ALSO CALCULATES THE CENTER OF GRAVITY OF TOE AIRCRAFT , 
THE CENTER OF GRAVITY SHIFT DURING THE TRAJECTORY MISSION AND THE 
FORWARD AND AFT CENTER OF GRAVITY LIMITS IMPOSED BY THE SIZE OF THE 

HORIZONTAL CONTROL SURFACE. , nc TUC 

IT ALSO DETERMINES THE PITCHING MOMENT CURVE SLOPE (DCM/DCL) OF THE 
AIRCRAFT FOR EACH PHASE IN THE MISSION. 

THE REFERENCES SITED THROUGHOUT THE PROGRAM ARE LISTED AT THE END 
OF THE LISTING. 


THE FILE ’STBCM . INC ' MUST BE IN THE SAME DIRECTORY AS THIS MODULE 
FOR IT TO WORK. THE FILE CONTAINS ALL THE COMMON BLOCK STATEMENTS 
NEEDED TO RUN THIS MODULE. THE 'INCLUDE’ STATEMENT IS USED TO 
INCLUDE THE FILE IN WITH THE MODULE DURING COMPILATION. S.M.S. 1/89 

THE FOLLOWING IS A LISTING OF THE STBCM. INC FILE 

*********** 

THIS CONTAINS THE COMMON BLOCK STATEMENTS NEEDED TO RUN THE STABILITY 
AND CONTROL MODULE ’STBLCON.FOR’ . THIS FILE IS C0NN j£TED WITH THE 
MODULE THROUGH THE ’INCLUDE’ STATEMENT FOUND AT THE BEGINING OF EACH 
OF THE SUBROUTINES. 

COMMON/OVER/ICALC , NERROR , MODNO , IGEO ,KGPRNT, IGPLT , IF 

COMMON/STBCM/CLALFA , CLAFHT , DEDA , XqCBCN , XQCBWG , XRCBHT , XQCBVT , 

1 CBARCN , CBARWG , CBARHT , CBARVT , ROOTWG , TRWG , XLEWG , SPANWG , 

2 SWPWG , POOL , EN , BODL , XLEPOD (10) , WCAND , WFTOT , FUFRAC , WWING , WFS , 

3 WLC , WHT , WVT , WETANK , WFEXT , WBODY , WAIRC , WAPU , WSC , WENG , SCN SWC SHT 

4 SVT WCREW , WELT , WEP , WHDP , WINST , WNA , WPA , WCARGO , WAMMUN , WBOMB , WMISS , 

5 ARWC,DIHED,ZRTWG,ZRTHT ,ARVT ,SWPCN,TRCN,AMTO, 

6 ZRTCN,VOLB,BDMAX,FRN,FRAB, ALPHA, CL, THRUST, ARHT,TNT(12) ,WFT(12) , 

7 WFTO WTOT ARCN , CMWB , WGTO , WFEO , WBB2 , WAF , WE , WPL , WPS , WTSUM , WFUEL , 

8 TWTO , MACH , ALT , SWETWG , SWETHT , SWETCN , SPANHT , SPANCN , 

9 STARTM(12) , DUMMY, IAO 

COMMON/LOC ALS/ETAHT , XFC AN , XFHT , XFNA , XFVT , XFWO , XFAFLLXFENG 

1 XFELT XFEP XFINST,XFHDP,XFPA,XFPAYL,XFFFUS,XFCREW,WPAYL,XFFUR, 

2 XFLCFRT , WFLGFRT , XFLGR , XFFS , XFBOMB , XFSC , XFPTV , XFLIFTF , WPT V WFUR , 

2 XEF1 , XFAMMUN , XFBB2 , XFMISS , XFC ARGO , XACWG , ETACN , ODEDA , XAWB , CLACN , 

3 WZET WZEF WEFF , WEFW , XEF2 , WEFWT , WZENC , SUMMO , SUMWT , XLE , LWB , LTTWB , 

3 IDBPRT , XCG (13) , DCMDCL (13) , SM (13) , XCGAFT (13) , WLCFRT , JjLCJ * » > 

4 WFP , WZB , WZM , WZ A , IBS , IMS , IAS , NPHASE , ICCPRT , WLFTF , WFFUS , LHT LCN , 

5 CELV CBALHT , TZCHT , FNEUTPT (13) , XFBOD , XFFWG , XFEX , XLETIP , IETANK , 

6 ZMTOR, AWTOT, ZMU.XTl ,XT2, ZFMG , ZMG , ZFT1 , ZTl , ZFT2, ZT2, ZFD ,ZD ,T1 , 

7 TSPLIT,CAMAT1,GAMAT2,IHT,ICN,ELVDEF,VR0T,ALPHACN,TAUJ2 

8 SHT 1 , SCN1 , ELVDMAX , IWG , XCGFRWD (13) ,QT0R,CGBAR(13) ,ALFTO,CLTO, 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


0 CDTO , CMWBTO , ALFMAX , CDMAX , CMWBMAX , ALFO CLO ' jj° ’ ScRtWVECT 

9 T1RAD , T2RAD , XLECG , XLET1 , X^ET2 , CLMAX , CLAHTGW , CLACNGRD - 

1 SPANCN1 , ZRTCN1 , ARCN1 , TRCN1 , ALFHT , ALFCN , DCMTMUL , TWZTO 

COMMON/XCGS/XBOD , XWING , XCAN , XHT XVT XNA XPIV XAIRC XAPU , 

1 XELT XEP XIN5T , XHDP , XSC , XPA , XENG ( 10) , XFS , XLIFTF , XFUR , 

2 XPAYL , XCREW , XCARCO , XAMMUN , XBB2 , XBOMB , XMIS5 , XETANK , 

3 XFUSF , XWGFUEL , XLGFRT , XLGR , LFLAG 

dcai YArwr DfMDCLF DCMDCLT . DCMDCLW , DCMCLCN , KF , CSC , TC , MBOD , 

1 MCAN , MWG , MFS , MLG , MHT , MVT , MET ANK , ^EXT , MAIRC , 

, 1 WR I Hi L H2 LH3 LH4 LH5 , LH6 f LH7 , LH8 ,LH9,LC1 , LC2 r LC3 ,LC4 f LLb, 

2 LC6 LC7 LC8 LHT LCN ', ZMt6r , ZMU , LTTWB , IWG , IHT , ICN , LEFF , MACH 

2 MCREW MELT, MEP.MHDP, MINST, MNA, MPA, MCARGO.MAMMUN MBOMB,MMISS, 

3 MPAYL , MPIV , MLIFTF ,K11,K12,K13,K14, K21 , K22 , K23 , K24 


INCLUDE ’STBCM . INC ’ 


IF (ICALC.GT.l .) GO TO 10 
C READ IN THE INPUT PARAMETERS. 

CALL STABINPT 

C 

10 CONTINUE 

IF (ICALC.NE.2.) GO TO 20 
C DO ALL THE NECESSARY CALCULATIONS. 

CALL STABCALC1 
C 

20 CONTINUE 

C IF (ICALC.NE.3.) GO TO 30 

C DO A FINAL PRINTOUT OF THE DATA. 
CALL STABOUT 

C 

30 CONTINUE 

C 

RETURN 

END 


* 










SUBROUTINE STABINPT 

THIS SUBROUTINE READS IN THE INP"T FROM THE INn" fREATEP BY 

tup 1 icpR TN NAMELIST 'STABIN ? IT ALSO INITILIZEo ALL VALUES SO 
THE USER DOES NOT NEED TO SPECIFY Y ALL VALUES, ONLY THE ONES DESI E 

INCLUDE ’STBCM. INC’ 

NAMELIST /STABIN/ETAHT XFC AN , XFHT , XFNA , XFVT , XFWG , XFAPU , XEBB2 , 

1 xFPIV XFELT XFEP XFINST,XFHDP,XFPA,XFFFUS , XFBOD , XFFWG , XFCARGO , 

2 XFLGFRT , WFLGFRT , XFLGR , XFFUR . XFMTSS , XFAMMUN . XFPAYL , XFEX , XFFc- , 
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C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


2 XFSC , XFBOMB , XFENG , XFLIFTF , XFCREW , IDBPRT , CELV , CBALHT , TZCHT , 

3 ETACN , AWTOT , ZMU , XT1 , XT2 , ZFMG, ZFT 1 , ZFT2 , ZFO , TSPLIT , GAM ATI , 

4 GAMAT2 , IHT, ICN , ELVDEF , ELVDMAX , VROT , SCN1 , XLECN , 

5 SPANCN1 , ZRTCN1 , TRCN1 , IETANK , IBS , IMS , IAS , NPHASE , ICGPRT , IVECT 


THE FOLLOWING IS A LIST OF THE INPUTS AND WHAT THEY STAND FOR... 

ETAHT - AFT HORIZONTAL TAIL EFFICIENCY 

ETACN - FORWARD CANCARD EFFICIENCY 

XFBOD - FUSSELAGE BODY C.G. POSITION IN PERCENT BODL 

XFCAN - FORWARD CANARD C.G. POSITION IN PERCENT MAC 

XFHT - AFT HORIZONTAL TAIL C.G. POSITION IN PERCENT MAC 

XFNA - NACELL(S) C.G. POSITION IN PERCENT BODL 

XFVT - VERTICAL TAIL C.G. POSITION IN PERCENT MAC 

XFWG - WING C.G. POSITION IN PERCENT MAC 

XFLGFRT - NOSE GEAR C.G. POSITION IN PERCENT BODL 

WFLGFRT - PERCENT NOSE GEAR OF TOTAL LANDING GEAR WEIGHT 

XFLGR - MAIN LANDING GEAR C.G. POSITION IN PERCENT BODL 

XFPIV - CONTROL PIVOTS C.G. POSITION IN PERCENT BODL 

XFFUR - FURNISHINGS C.G. POSITION IN PERCENT BODL 

XFAPU - APU C.G. POSITION IN PERCENT BODL 

XFELT - AVIONICS C.G. POSITION IN PERCENT BODL 

XFEP - ELECTRICAL SYSTEMS C.G. POSITION IN PERCENT BODL 

XFINST - INSTRUMENTS C.G. POSITION IN PERCENT BODL 

XFFS - FUEL SYSTEM C.G. POSITION IN PERCENT BODL 

XFHDP - HYDRAULIC AND PNEUMATIC C.G. POSITIONS IN PERCENT BODL 

XFSC - SURFACE CONTROLS C.G. POSITION IN PERCENT BODL 

XFPA - PASSENGER ACCOMADATIONS C.G. POSITION IN PERCENT BODL 

XFPAYL - PAYLOAD C.G. POSITION IN PERCENT BODL 

XFCARGO - CARGO C.G. POSITION IN PERCENT BODL 

XFFWG - FUEL IN WING C.G. POSITION IN PERCENT BODL 

XFFFUS - FUEL IN FUSELAGE C.G. POSITION IN PERCENT BODL 

XFEX - EXTERNAL TANKS C.G. POSITION IN PERCENT BODL 

XFCREW - CREW C.G. POSITION IN PERCENT BODL 

XFENG - ENGINE (S) C.G. POSITION IN PERCENT PODLength 

XFLIFTF - LIFT FAN C.G. POSITION IN PERCENT BODL 

XFBOMB - BOMBS C.G. POSITION IN PERCENT BODL 

XFBB2 - ADVANCED WEAPONS SYSTEMS C.G. POSITION IN PERCENT BODL 

XFMISS - MISSILES C.G. POSITION IN PERCENT BODL 

XFAMMUN - AMMUNITION C.G. POSITION IN PERCENT BODL 

IDBPRT - DEBUG PRINT FLAG, CREATES MORE OUTPUT 

CELV - PERCENT OF STABILIZER WHICH IS MOVABLE ELEVATOR 

CBALHT - PERCENT OF ELEVATOR WHICH IS FORWARD OF THE HINGE LINE 

TZCHT - THICKNESS TO CHORD RATIO OF THE ELEVATOR 

AVYTOT - ANGLE OF ATTACK OF THE WING AT TAKE OFF ROTATION 

ZMU - FRICTION COEFFICIENT OF RUNWAY ON TTRES 

XU - HORIZONTAL FOSITTON OF FORWARD THRUST VECTOR PERCENT BODL 

XT2 - HORIZONTAL POSITION OF AFT THRUST VECTOR IN PERCENT BODL 

ZFMG - VERT DIST. FROM FUSS. CENTERLINE TO MAIN GEAR, % BDMAX 

ZFTl - VERT DIST. FROM FUSS. CENTERLINE TO FORE THRUST, 3BDMAX 

ZFT2 - VERT DIST. FROM FUSS CENTERLINE TO AFT THRUST, JtBDMAX 

ZFD - VERT. DIST. FROM FUSS. CENTERLINE TO DRAG VECTOR , /5BDMAX 

TSPLIT - FORWARD THRUST TO TOTAL THRUST RATIO IN PERCENT 

CAMAT1 - FORWARD THRUST VECTOR ANGLE, DEGREES 
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C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

r 

c 

c 

c 


GAMAT2 

SCN1 

XLECN 

SPANCN1 

ZRTCN1 

TRCN1 

IHT 

ICN 

ELVDEF 

ELVDMAX 

VROT 

IETANK 


IBS 

IMS 

IAS 

NPHASE 

ICGPRT 

IVECT 




AFT THRUST VECTOR ANGLE, DEGREES 

CANARD AREA USED WHEN A.C. HAS A AFT TAIL ALSO SQ_ FT. 
CANARD LEADING EDGE POSITION USED WITH AFT TAIL, FEET 
CANARD SPAN USED WITH AFT TAIL, FEET 
VERTICAL CANARD HEIGHT USED WITH AFT TAIL, FEET 
TAPER RATIO OF CANARD USED WITH AFT TAIL 
AFT TAIL INCEDENCE ANGLE, DEGREES 
FORWARD CANARD INCEDENCE ANGLE, DEGREES 
ELEVATOR ANGLE NEEDED FOR TAKE OFF ROTATION, DEGREES 
MAXIMUM ELEVATOR DEFLECTION FOR LANDING, DEGREES 
MACH NUMBER AT TAKEOFF ROTATION 
EXTERNAL TANK FLAG, USED TO SPECIFY POSITIONS OF 
ANY EXTERNAL TANKS. 0=TANK ON FUSS. CENTERLINE, l=TWO 
FUSS SIDE MOUNTED TANKS, 2=TWO MAC MOUNTED TANKS, 

3=TW0 WING-TIP TANKS, 4=ONE CENTERLINE TANK AND TWO 
MAC MOUNTED TANKS, 5=CENTERLINE TANK AND WING-TIP TANKS, 
6=TW0 MAC TANKS AND TWO WING-TIF TANKS. 

INDICATOR TO REMOVE BOMB WEIGHT ON A SPECIFIED MISSON PHASE 
INDICATOR TO REMOVE MISSILE WEIGHT ON SPECIFIED PHASE 
INDICATOR TO REMOVE AMMUNITION WEIGHT ON SPECIFIED PHASE 
NUMBER OF PHASES IN MISSION 

PRINT FLAG TO PRINT C.G. INFORMATION FOR EACH PHASE 
VECTORED THRUST INDICATOR, 0= SOLVE FOR NO VECTORED 
THRUST PARAMETERS, 1=S0LVE FOR FORWARD THRUST VECTOR 
POSTION, 2=S0LVE FOR THRUST SPLIT, 3=S0LVE FOR AFT 
THRUST VECTOR ANGLE. 
******»«**********'»***<"***********'**** , *' , '****** , ****'* , '"*** l * , ' , ** ,, '"‘' 


INITIALLZE THE STABIN INPUT VALUES 


ETAHT =0.90 
ETACN =0.90 
XFBOD = .50 
XFCAN = .20 
XFHT = .20 
XFNA = .40 
XFVT = .20 
XFWG = .15 
XFLGFRT = .10 
WFLGFRT = .3 
XFLGR = .55 
XFPIV = .65 
XFFUR = .50 
XFAPU = .95 
XFELT = .05 
XFEP = .50 
XFINST = .15 
XFFS = .60 
XFHDP = .50 
XFSC = .70 
XFPA = .55 
XFPAYL = .50 
XFCARGO = .50 
XFFWG = .50 



n n n n 


74 


XFFFUS = .50 
XFEX = .50 
XFCREW = .20 
XFENG = .50 
XFLIFTF = .50 
XFBOMB = .50 
XFBB2 = .50 
XFMISS = .47 
XFAMMUN = .30 
IDBPRT = 0 
CELV =0.30 
CBALHT =0.0 
TZCHT =0.10 
AWTOT =2.0 
ZMU = .02 
XT1 = .30 
XT2 = .60 
ZFMG = 1 .0 
ZFT 1 = 0.0 
ZFT2 = 0.0 
ZFD = 0.0 
TSPLIT =0.0 
GAMAT1 = 0.0 
GAMAT2 =0.0 
SCN1 = 0.0 
XLECN =0.0 
SPANCN1 =0.0 
ZRTCN1 =0.0 
TRCNl =0.0 
IHT = 0.0 
ICN = 0.0 
ELVDEF = -10.0 
ELVDMAX = -15.0 
VROT =0.3 
IETANK = 0 
IBS = 0 
IMS = 0 
IAS = 0 
NPHASE = 1 
ICGPRT = 0 
IVECT = 0 

READ IN THE VALUES FROM THE NAMELIST 

IF (ICALC .CT. 1 .) GO TO 2°0 
READ (5, 40) TITLE 
READ(5,STABIN) 

OUTPUT THE INPUTS TO THE OUTPUT FILE SO THE USER CAN DOUBLE 

CHECK THAT THE VALUES INPUT ARE THOSE DESIRED 

WRITE (6, 300) 

WRITE (6 , 310) ETAHT , XFBOD , XFCAN 

WRITE (6,315) ETACN , WFLGFRT 



WRITE (6 , 320) XFHT , XFNA , XFVT 
WRITE (6 , 330) XFWG , XFLGFRT , XFLGR 
WRITE (6 , 340) XFPIV , XFFUR , XFAPU 
WRITE (6 , 350) XFELT , XFEP , XFINST 
WRITE (6 , 360) XFFS , XFHDP , XF5C 
WRITE (6 , 370) XFPA , XFPAYL , XFFWG 
WRITE (6 , 380) XFFFUS , XFEX , XFCREW 
WRITE (6 , 390) XFENG , XFLIFTF 
WRITE (6 , 400) XFAMMUN , XFBOMB , XFMISS 
WRITE (6 , 405) XFCARGO , XFBB2 
WRITE (6,410) CELV , CBALHT , TZCHT 
WRITE (6 , 4 1 1 ) AWTOT , ZMU 
WRITE(6,412)XT1,XT2,ZFMG 
WRITE (6 , 413)ZFT1 , ZFT2 ,ZFD 
WRITE (6 , 4 1 4) TSPLIT , CAMAT 1 , GAMAT 2 
WRITE(6,415)IHT,ICN,ELVDEF 
WRITE (6,416) VROT , SPANCN1 , ZRTCN1 
WRITE (6 , 417) SCN1 ,XLECN,ELVDMAX 
WRITE(6 , 418)TRCN1 , IVECT 
WRITE (6 , 420) IETANK , IBS , IMS , IAS 
WRITE (6 ,430) NPHASE , IDBPRT , ICGPRT 
290 CONTINUE 


40 FORMAT (20A4) „ 

300 FORMAT (// , 1H1 , 15X,28HSTABILITY AN ^ CONTROL J^S,//) F 0 5> 

310 FORMAT (10X,7HETAHT= ,F10.5,5X,8HXFBOD= ,F10.5,5X,7HXFCAN= ,F10.5) 

- , F10 . 5 , 5X , 8HWFLGFRT= , F10 . 5) 

F10 . 5 , 5X , 8HXFNA= ,F10.5,5X,7HXFVT= ,F10.5) 

F10 . 5 , 5X , 8HXFLGFRT= , F10 . 5 , 5X , 7HXFLGR= , F10 . 5) 
*F10 . 5 , 5X , 8HXFFUR= ,F10.5,5X,7HXFAPU= ,F10.5) 
[ FlO . 5 , 5X , 8HXFEP= 

, F10 . 5 , 5X , 8HXFHDP= 

, v _.. f , F10 . 5 , 5X , 8HXFPAYL= 

380 FORMAT (10X j 7HXFFFUS= , F10 . 5 , 5X , 8HXFEX= . 

390 FORMAT (10X, 7HXFENG= ,F10.5,5X,8HXFLIFTF=,F10.5) 

400 FORMAT (10X , 7HXFAMUN= , F10 . 5 , 5X , 8HXFB0MB= , F10 . 5 , 5X , 7HXFMISS= , FlO . 5) 

405 FORMAT (10X,7HXFCARG=, FlO. 5, 5X,8HXFBB2 = C1rt K ' 

410 FORMAT (10X , 7HCELV= ,F10.5,5X,8HCBALHT= 


315 FORMAT (10X,7HETACN= 
320 FORMAT (10X,7HXFHT= 
330 FORMAT (10X,7HXFWG= 
340 FORMAT (10X , 7HXFPIV= 
350 FORMAT (10X,7HXFELT= 
360 FORMAT (10X,7HXFFS= 
370 FORMAT (10X,7HXFPA= 


111 FORMAT (10X r 8HAWTOT = , FlO . 5 , 5X, 7HZMU = ,F10.5) 

_ x ~ ^ v/ ftiivTn _ cm c 


, FlO . 5 , 5X , 7HXFINST= , FlO 5) 
, F10.5,5X,7HXFSC= ,F10.5) 
, FlO . 5 , 5X , 7HXFFWG= ,F10.5) 
, FlO . 5 , 5X , 7HXFCREW= , FlO . 5) 


,F10.5) 

,F10.5,5X,7HTZCHT= ,F10.5) 


412 FORMAT (10X,7HXT1 = ,F10.5,5X,8HXT2 = 

413 FORMAT (10X,7HZFT1 = ,F10.5,5X,8HZFT2 = 

4 1 4 FORMAT ( 1 OX , 7HTSPLIT= , F 1 0 . 5 , 5X , 8HG AMAT 1 = 

415 FORMAT (10X,7HIHT = ,F10.5,5X,8HICN = 

416 FORMAT (10X,7HVROT = ,F10.5,5X,8HSPANCN1 


, FlO . 5 , 5X , 7HZFMG = ,F10.5) 
, F10.5,5X,7HZFD = ,F10.5) 
, FlO . 5 , 5X , 7HGAMAT2= , F 10 . 5) 
,F10.5,5X,7HELVOEF=,F1.0.5) 
FlO 5 ,5X, 7HZRTCN1 = , FI 0 . 5) 


17 FORMAT (10X,7HSCN1 = ,F10.S,5X,8HXLFCN = , F 1 0 5 , 5X , 7HFI_VI)MX= , O 

18 FORMAT (10X,7HTRCN1 = ,F10 5 , 5X, 8HTVECT = ,13) 

20 FORMAT (10X, 7HIETANK= , 13 , 5X, 5HIB3= ,I3,5X,5 H IMo= 13 5X,5HIAS= ,13) 

30 FORMAT (10X, 7HNPHASE= , 13 , 5X, 8HIDBPRT= ,I3,5X,8HICGPRT_ ,13) 


RETURN 

END 
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SUBROUTINE STABCALC1 
THIS SUBROUTINE DETERMINES THE 

SagSeaSaSSsaB* 

EsF sbf 

PHASES AND THE WEIGHT LOSSES DUE TO THE USE OF WEAPONS. 

INCLUDE ’STBCM INC ’ 




C 

c 


I FLAG = 0 
CGFLAG = 0 
LCFLAG = 0 
112 CONTINUE 
T = 1 

IPHASE = 0 
XCGFRWD (1) = 0 

ciSSTthf acsynt calculated weights in local variables and 

sns ™1 EUEL^WEIGHT TO THE MAXIMUM VALUE FOR TAKE-OFF 
CONDITIONS. 

WFP = WFTOT 
WZB = WBOMB 
WZM = WMISS 
WZA = WAMMUN 
WZET = WETANK 
WZEF = WFEXT 

SIZE THE TAIL FOR TAKE-OFF ROTATION IF IT HAS NOT ALREADY BEEN DONE 
IF (I FLAG EQ. 0 AND. CGFLAG .EG. 0) THEN 

CALL TAILSIZE (I, IPHASE) 

END IF 
I FLAG = 0 
CGFLAG = 0 

Tck ^SfW?5sS “-ot Boreas ™”ua 
$ ™e i ectIInaL ™nks if it IS the tanks are removed 

1=2 

IF (WFEXT . GT. 0.) THEN 
WFP = WFP - WFEXT 
IF (WFTO .CE. WFEXT) THEN 
WZET = 0. 

WZEF = 0. 

WFP = WFP - (WFTO-WFFXT) 


ELSE 


WZEF = WFEXT - WFTO 



0 0 0 0 
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c 


ENDIF 

ELSE 

WFP = WFP - WFTO 
END IF 


CALL STABCALC2 (I , IPHASE , CGFLAG) 
IF (CGFLAG .EH. 1)THEN 
CO TO 112 
END IF 


REMOVE THE FUEL WEIGHT USED IN EACH PHASE, AND REMOVE ANY WEAPONS 
USED DURING EACH PHASE. ALSO REMOVE EXTERNAL TANKS WHEN THE FUEL 
USAGE BECOMES GREATER THEN THE SIZE OF THE EXTERNAL TANKS. 


DO 100 IPHASE = 1 ,NPHASE 
IF (WZEF .GT. 0.) THEN 

IF (WFT (IPHASE) .CE. WZEF) THEN 

WFP = WF p _ (WFT (IPHASE) -WZEF) 
WZET = 0. 

WZEF = 0. 

ELSE 

WZEF = WZEF - WFT (IPHASE) 

ENDIF 

ELSE 

WFP = WFP - WFT (IPHASE) 

ENDIF 

IF (IPHASE. GE. IBS) WZB = 0. 

IF (IPHASE. GE.IWS) WZM = 0. 

IF (IPHASE . GE. IAS) WZA = 0. 

C 

I = IPHASE ♦ 2 


CALL STABCALC2 (I, IPHASE, CGFLAG) 
IF (CGFLAG .EQ. 1)THEN 
GO TO 112 
ENDIF 


C CALCULATE THE FORWARD C.G. LIMIT FOR THE TAIL SIZE DETERMINED EARLIER. 
C 

CALL XCGFORWARD(I) 


THESE SET OF INSTRUCTIONS DETERMINE IF THE CENTER OF GRAVITY I' TKirD(rA . r . 
BEHIND THE FORWARD C.G. LIMIT DETERMINED, AND IF IT IS NOT, THEN TT INCRFA r 
THE SIZE OF THE HORIZONTAL CONTROL SURFACE AND RECALCULATE..' THE 
FORWARD C.G. LIMIT. 


IF (XCGFRWD (I) .GT. XCG(I.))THEN 
WRITE (6, 300) XCG (I) ,XCCFRWD(I) 
IF (SCN .EG. 0.0) THEN 
SHTl = SHT1 * 5.0 
I FLAG = 1 
ELSE 

SCN1 = SCN1 ♦ 5.0 
I FLAG = 1 
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ENDIF 
GO TO 112 
ENDIF 
C 

100 CONTINUE 

c THIS SET OF INSTRUCTIONS ACTIVATES THE VECTORED THRUST ANALYSIS 

C IF IVECT IS NOT SET TO 0. 

C 

IF (IVECT .EQ. 0) GO TO 200 
J = NPHASE + 2 
CALL VECTHRUST (I) 

200 CONTINUE 


C 


300 FORMAT (/ , 10X , 35HWARNINC ! ! C.C. IS FORWARD OF LIMI ^M 5X > 

139H CONTROL SURFACE INCREASED BY 5 SQ. FT . , / , 15X , 6HXCG = ,F12.6, 
25X , 10HXCGFRWD = ,F12.6) 


RETURN 

END 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

r 

r 
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SUBROUTINE TAILSIZE(I , IPHASE) 

THIS SUBROUTINE IS USED TO CALL THOSE ROUTINES NEEDED TO SIZE THE 
HORIZONTAL TAIL AT TAKEOFF. 

INCLUDE ’STBCM . INC ’ 

DATATRANS SUBROUTINE TRANSFERS OVER THE NECESSARY AERODYNAMICS 
PARAMETERS . 

CALL DATATRANS 

THIS SUBROUTINE DOES ALL THE SECONDARY AERODYNAMICS CALCULATIONS. 
CALL STABCALC3 (I , IPHASE) 

THIS SUBROUTINE DETERMINES THE TAIL SIZE NEEDED FOR ROTATION. 

CALL TOROTATE (I, IPHASE) 

RETURN 

END 


C 

C 

c 

c 


SUBROUTINE STABCALC2 (I , IPHASE , CGFLAG) 

THIS SUBROUTINE IS USED TO CALCULATE ANY PARAMETERS NEEDED ArcvMT 
BY THE SUBROUTINES WHICH ARE NOT CALCULATED BY THE REST OF ACSYNT 

INCLUDE ’STBCM. INC’ 
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THIS SERIES OF INSTRUCTIONS SETS THE MACH NUMBER EITHER TO THE 
TAKE OFF MACH NUMBER OR THE MACH NUMBER AT THE START OF A MISSION 
PHASE. IT THEN CALLS AERODYNAMICS TO GET THE NEW CLALPHA OF THE 
WING AND THEN CALLS THE AERODYNAMIC CENTER SUBROUTINE. 

IF (IPHASE .EQ. 0) THEN 
IF (AMTO . LE. 0.0) THEN 
AMTO = VROT 
ENDIF 

MACH = AMTO 
ELSE 

MACH = STARTM (IPHASE) 

ENDIF 
ICALC = 2 
IAO = 1 

CALL STBDT (ICALC , NERR , IGEO , KGPRNT , IGPLT) 

2CLALFA = CLALFA 
ZMACH = MACH 

NOW CALCULATE THE AERODYNAMIC CENTER INCORPERATING MACH EFFECTS, 
SWEEP ANGLE EFFECTS, AND TAPER EFFECTS. 

CALL XACCALC (ZMACH, ZCLALFA) 

NOW CALCULATE THE CENTER OF GRAVITY OF THE AIRCRAFT. 

CALL XCGCALC (I) 

200 CONTINUE 

THE FOLLOWING SERIES OF COMMANDS CALCULATE THE PITCHING MOMENT 
CURVE SLOPE FOR THE AIRCRAFT. THESE CALCULATIONS ARE BASED ON 
EQUATIONS TAKEN FROM REFERENCE 4, CHAPTER 5. THEY INCLUDE THE 
TAIL, WING, CANARD, FUSELAGE, AND JET(S) . 

.♦TAIL CONTRIBUTION TO AIRCRAFT STABILITY 

CSC=CLALFA.SWG*CBARWG 

DCMDCL1=-CLAFHT*SHT 1 * (XQCBHT-XCG (I) ) .ETAHT . (1 . -DEDA) /CSC 
THIS ACCOUNTS FOR STICK FREE EFFECTS, DETERMINED IN STABCALC3 
DCMDCLT=DCMDCL1 .DCMTMUL 

..CALCULATE THE PITCHING MOMENT SLOPE CONTRIBUTION DUE TO 
THE CANARD NOTE THAT THERE ARE TWO SETS OF COMMANDS FOR 
THE CANARD, ONE FOR A CANARD ALONE, AND ONF FOR A CANARD 
AND AFT TAIL CONFIGURATION 

TF (XLECN .EQ. 0.0) THEN 

VOLCN = (SCN1/SWG).((XCG(I)-X0CBCN)/CBARWG1 

ELSE 

VOLCN = (SCN1/SWC).((XCG(I)-XLECN)/CBARW0) 

ENDIF 
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DCMCLCN=(CLACN/CLALFA)*VOLCN*ETACN*ODEDA 

COMPUTE THE STICK FREE DCMDCL OF THE CANARD IF THERE IS 
NO H0R0Z0NTAL STABALIZER, OTHERWISE THE H.T. IS ASSUMED 
TO BE FREE AND THE CANARD IS USED ONLY FOR TRIM. 

IF (SCN CT. 0. .AND. SHT .LE. 0.) THEN 
DCMCMUL = DCMTMUL 
DCMCLCN = DCMCLCN * DCMCMUL 
ELSE 

DCMCMUL = 1. 

END IF 

* *WING CONTRIBUTION TO AIRCRAFT STABILITY 
DCMDCLW= ( (XCG(I) -XLE) /CBARWG) - ( (XACWG-XLE) /CBARWG) 

..FUSALAGE CONTRIBUTION TO STABILITY 

THIS USES THE EASY METHOD FROM REFERENCE 4 WHICH CURVE-FITS 
THE GRAPH OF KF (FIGURE 5-16) TO DETERMINE THE FUSSELAGE EFFECTS. 

KF= ( (XLEWG+R00TWG/4 . ) /BODL) *100 . 

KF= 0023+ .0002062.KF- .000005762*KF»*2+.0000002499»KF»*3 
CLARAD = CLALFA * 57.29578 

DCMDCLF=(KF* (BDMAX**2) *BODL) / ( SWG * CBARWG * CL ARAD) 

CALCULATE THE AERODYNAMIC CENTER FOR THE WING AND FUSSELAGE BY 
ADDING THE EFFECTS OF THE FUSELAGE TO THE AERODYNAMIC CENTER 
OF THE WING. 

XAWB = ((XACWG-XLE) /CBARWG) * DCMDCLF 

.♦STABILITY EFFECTS DUE TO THE JET(S) 

1; EFFECTS DUE TO THRUST 

TNI = TNT (IPHASE) *TSPLIT 

TN2 = TNT (IPHASE) - TNI 

DDT1 = (TNI * ZTl) / (SUM WT * CBARWG) 

DDT2 = (TN2 * ZT2) / (SUMWT * CBARWG) 

DDCMTH = DDTl ♦ DDT2 

2; EFFECTS DUE TO NORMAL FORCE (THERE ARE NOT ANY YET.) 

DDCMNF = 0. 

TOTAL JET EFFECTS 
DCMDCL J = DDCMTH + DDCMNF 

..OVERALL MOMENT-CURVF SLOPE 

THIS IS JUST A SUMMATION OF ALL THE PARTS. 

DCMDCL (!)=DCMDCl.W+PCMCLCN+RCMDCLF+DCMDCLT*ncMDCLJ 

..AIRCRAFT STICK-FREE HP URAL POINT 
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c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 


THIS USES 
MODIFIED 
FUSELAGE 


AN EQUATION FROM REFERENCE 2 (CHAPTER 5-EQN. 5.145) 
TO INCLUDE A CANARD. ASSUMES THRUST EFFECTS AND 
EFFECTS DO NOT CHANGE WITH ANGLE OF ATTACK. 


CLHR = (CLAFHT /CLALFA) * (SHT 1 /SWG) *ETAHT 
CLCR = (CLACN/CLALFA) * (SCN1/SWG) *ETACN 

FNEUTPT (I) = (XAwScLHR* (XRCBHT/CBARWG) ♦ (1 . -DEDA) *DCMTMUL) / (1 . 
S +CLHR* (1 . -DEDA) *DCMTMUL) 

ELSE 

FNEUTPT (I) = (XAWB+CLCR* (XQCBCN/CBARWG) .ODEDA.DCMCMUL) / (1 • ♦ 

S CLCR* ODEDA.DCMCMUL) 

END IF 


»»CG AFT LIMIT 

THE AFT C.G. LIMIT IS THE NEITTRAL POINT OF THE AIRCRAFT. 
XCGAFT (I + l) = XLE ♦ (FNELTTPT (I) * CBARWG) 


..STATIC MARGIN 

THE STATIC MARGIN IS SIMPLY THE DIFFERENCE BETWEEN THE 
NEUTRAL POINT AND THE CENTER OF GRAVITY. 


SM(I)= ( (XCGAFT (1*1) - XCG (I)) /CBARWG) *100. 


CHECK TO SEE IF THE C.G. FOR EACH PHASE IS AFT OF THE AFT C G 
LIMIT IF IT IS, THEN THE TAIL SIZE IS INCREASED, AND THE 
ENTIRE SUBROUTINE RERUN 


IF (XCGAFT (I+l) .LT. XCG(I))THEN 
WRITE (6, 300) XCG (I) ,XCGAFT(I*1) 
IF (SCN .EQ. 0.0) THEN 
SHT 1 = SHT1 ♦ 5.0 
CGFLAC=1 
RETURN 
ELSE 

SCN1 = SCN1 * 5.0 
CGFLAG=1 
RETURN 
END IF 
END IF 


IF (IDBPRT.EQ.O) GO TO ac,o 


THIS IS THE WRITE STATEMENT? TO BE USED WHEN COMPLETE OEBUGING 
INFORMATION IS DESIRED. THIS IS STARTED BY SETTING IDBPRT TO 1 


IF (IPHASE .LE. 0)THEN 
WRTTE(6,59) 

ELSE 

WRITE (6, 60) IPHASE 
ENDIF 
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WRITE (6, 100) 
WRITE (6, 105) 
WRITE(6, 110) 
WRITE (6, 120) 
WRITE (6, 130) 
WRITE(6,140) 
WRITE (6, 410) 
WRITE(6,420) 
WRITE (6 , 425) 
WRITE (6 , 430) 
WRITE (6, 440) 
WRITE (6, 450) 


XCC(I) 

(XAWB*CBARWG) +XLE 
DCMDCL (I) 

FNEUTPT (I) 

SM(I) 

XCGAFT(I*1) 

SUMM0,5UMWT 

DCMDCLT,DCMDCLW 

DCMDCLF,DCMCLCN 

DCMDCLJ 

LWB, LHT , LCN 

SHT1 r SCNl 


400 CONTINUE 

59 FORMAT (5X,19HDATA AFTER TAKE-OFF) 

60 FORMAT (5X,2OHMISSI0N PHASE NUMBER, 15) 

100 FORMAT (10X,6HXCG = F15.7) 

105 FORMAT (10X,6HXAC = F15.7) 

110 FORMAT (10X,9HDCMDCL = F15.7) 

120 FORMAT (10X , 9HNEUTPT = F15.7) 

130 FORMAT (10X, 16HSTATIC MARGIN = F15.7) 

140 FORMAT (10X,9HAFT CG = F15.7) 

410 FORMAT (5X,8HSUMM0 = , FI 5 . 6 , 5X , 8HSUMWT = ,F15.6) 

420 FORMAT (5X, 10HDCMDCLT = , F15 8 , 5X, 10HDCMDCLW = ,F15.8) 

425 FORMAT (5X, 10HDCMDCLF = ,F15.8,5X,10HDCMDCLCN =,F15.8) 

430 FORMAT (5X, 10HDCMDCLJ = ,F15.8) 

440 FORMAT (5X,6HLWB = , F15 . 8 , 5X, 6HLHT = ,F15.8,5X,6HLCN = ,F15.8) 

450 FORMAT (5X,6HSHT = , F15 . 8 , 5X , 6HSCN = ,F15.8) 

300 FORMAT (10X, 31HWARNINC ! ! C.G. IS AFT OF LIMITS,/, 15X, 

144H CONTROL SURFACE SIZE INCREASED BY 5 SQ. FT. ,/,15X,6HXCG = , 
2F15 . 7 , 5X , 8HAFTCG = ,F15.7) 

RETURN 

END 


** * 


SUBROUTINE STABCALC3 (I , IPHASE) 

THIS SUBROUTINE DETERMINES THE AERODYNAMIC CENTER AND C.G. DURING 
TAKEOFF, IT DETERMINES THE STICK FREE EFFECTS, IT DETERMINES 
THE UPWASH EFFECTS ON THE CANARD, AND THE LIFT-CURVE SLOPE 
OF THE CANARD. 

INCLUDE ’STBCM. INC ’ 

DETERMINE THE AERODYNAMIC CENTER OF THE wtng, 

USING THE GEOMETRIC SHAPE AND THE 2-D AND 3-0 
LIFT COEFFICIENTS. THIS USES THE INPUT OF THE TAKEOFF MACH 
NUMBER . 

IF (AMTO .LE. 0.0) THEN 
AMTO = VROT 
ENDIF 

C NEED TO TRANSFER LIFT-CURVF SLOPE DATA FROM AERODYNAMICS MODULE 
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C TO DETERMINE THE AERODYNAMIC CENTER AT TAKE OFF MACH NUMBER. 

MACH = AMTO 
ICALC = 2 
IAO = 1 

CALL 5TBDT (ICALC , NERR , IGEO , KGPRNT , IGPLT) 

ZCLALFA = CLALFA 
ZMACH = MACH 
C 

CALL XACCALC (ZMACH, ZCLALFA) 

...» DETERMINE THE CENTER-OF-GRAVITY FOR THE AIRCRAFT •*** 

CALL XCGCALC (I) 

.* DETERMINE THE STICK FREE MUTIPLIER FOR THE HORIZONTAL TAIL. 

THIS WAS DONE BY CURVE FITTING TO DETERMINE THE TAIL 
HINGE MOMENTS DUE TO ALPHA AND ELEVATOR DEFLECTION. 

THE METHOD USED IS FROM REFERENCE 3, CHAPTER 8. EQNS 8.42 AND 
8.46 THE GRAPHS CURVE FIT ARE FIGURES 8.11, AND 8.12. 

INPUTS INCLUDE THE ELEVATOR SPAN, BALLANCE SPAN, 

AS WELL AS THE AR AND t/c OF THE TAIL. 

IF (SCN GT. 0.0 .AND. SHT .LE. 0.0) THEN 
CZELV = CELV.CBARCN 
CZBALHT = CBALHT . CELV 
CHT = CBARCN 
ARTAIL = ARCN 
THT = TZCHT ♦ CBARCN 

ELSE 

CZELV = CELV* CBARHT 
CZBALHT = CBALHT . CELV 
CHT = CBARHT 
ARTAIL = ARHT 
THT = TZCHT * CBARHT 

END IF 

C NEED TO DETERMINE THE BALLANCE RATIO OF THE ELEVATOR. IF IT IS 

C GREATER THEN 0.38, THEN IT IS SET TO 0.38 TO FIT THE GRAPHS. 

CF = CZELV - CZBALHT 

IF ((CZBALHT/CF)**2 . GT. (THT/(CF*2.))**2) THEN 
BR = SQRT (((CZBALHT/CF) **2) - ((THT/(CF*2.))**2)) 

ELSE 
BR = 0. 

END IF 

IF (BR .GT. 0.38) BR=0.?« 

CEC = CZELV/CHT 
A1A = 1. /ARTAIL 
TC = TZCHT 
r 

C HERE WE ARE CURVE-FITTING THE TWO CRAPHS TO GET THE VALUES FOR 
C K1 AND K2, BY WHICH WE THEN GET BGJ AND BC2 . 

C THESE ARE GRAPHS ARE FROM REFERENCE 3, FIGURE? 8.11 AND C.12 
C 

Kll=. 1027+3. 5772*CEC-1 ,8662*CEC**2* . 174*CEC*-3- ,3587-CEC**4 

K12=l . 1692 - 2 . 5003*TC*23.4O44*TC* *2-202. 18°P*TC*-*3*546 ,32S4*TC**4 
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K13=l . *2 . 8032*BR-31 . 5474 *BR **2+72. 841 2*BR* ♦3-56.8704 *BR* *4 
K14=l . -2 . 2462*A1A*1 . 2067‘A1 A**2-l . 378* Al A* *3*1 ,1478‘A1A**4 
BG1=-.55*K11*K12*K13*K14 

K21= . 7504* . 8193*CEC- .0829*CEC**2‘ .2082*CEC**3- . 1449*CEC**4 
K22=l .♦.15*BR-12.4349‘BR**2+18. 7076*BR**3-12 . 8205*BR**4 
K23=l . -1 . 0762*A1A+ . 41303»AlA*»2+ . 41321*A1A**3- . 37067 *A1A**4 
K24=l .0748- .7949+TC+9. 1 283 »TC»* 2-105 . 7365‘TC“3+233 .0985*TC“4 
BG2=- . 89*K21 *K22*K23‘K24 

TAU IS DETERMINED FROM REFERENCE 3, CHAPTER 3, FIGURE 3.32. 

THIS IS KNOWN AS THE ELEVATOR EFFECTIVNESS FACTOR. 

TAU=4 . 1065*CEC-8 .9175*CEC**2*9. 7178*CEC“3-3 . 9307*CEC**4 

HERE WE CALCULATE THE ’FREE ELEVATOR FACTOR’ (DCMTMUL) , WHICH WILL 
BE MULTIPLIED TO THE ELEVATOR ANGLE OF ATTACK TO DETERMINE STICK 
FREE EFFECTS. 

DCMTMUL=1- (TAU * (BC1/BC2)) 

* ‘CANARD CONTRIBUTION TO AIRCRAFT STABILITY 

FIRST DETERMINE THE UPWASH ON THE CANARD DUE TO THE WING. 

THIS IS FOUND BY CURVE FITTING FIGURE 8.22c FROM REFERENCE 3. 

CO = ROOTWG 

IF (XLECN .EQ. 0.0) THEN 

XFOR = (XLEWG * ROOTWC/4.) - XQCBCN 
ELSE 

XFOR = XLECN 
ENDIF 

XC = XFOR/CO 
IF (XC .CE. 5 . )THEN 
ODEDA = 1. 

GO TO 190 
ENDIF 

IF (ARWG .LE. 8.) THEN 
XI = 6. 

X2 = 8. 

Y1 =2 .6538-3 . 3623+XC+2 . 9508*XC**2 . -1 .2856‘XC*»3. 

S * . 2706*XC“4 ,-.0219»XC**5. 

Y2=2 . 8419-3 ,6494‘XC*3 . 1799»XC“2 . -1 . 3833*XC--3 
S *.2916*XC“4.-. 0236 *XC* *5 

ODEDA= ( (X2-ARWG) * Y 1 - (XI -ARWG) *Y?_) / (X2-X1) 

GO TO 190 
ENDTF 

IF (ARWG .GT. 8. .AND. ARWG .LE. 10.) THEN 
XI = 8. 

X2 = 10. 

Y1 =2 . 8419-3 . 6494 *XC *3 1799*X r “2.-1 ,3833*XC“3 
S * . 29I6*XC**4 . - . 0236‘XC**5 . 

Y2=2 . 9405-3 . 7399»XC*3 . 1866-XC* *2 -i 3592*XC“? 

S “2821‘XC“4-.0226 + XC“S 
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ODEDA= ( (X2-ARWG) *Y1 - (XI -ARWG) * Y2) / (X2-X1 ) 

GO TO 190 
ENDIF 

IF (ARWG .GT. 10.) THEN 
XI = 10. 

X2 = 12. 

Yl=2 .9405-3 . 7399*X03 . 1866*XC**2 . -1 . 3592*XC«*3 . 

$ + . 2821 *XC* *4 . - . 0226*XC»*5 . 

Y2=3. 1183-4 .1231 X"- .' . 5627*XC**2 . -1 . 5362*XC**3 . 

S + . 3217*XC**4 . - . 9259*XC** 5 . 

ODEDA= ( (X2-ARWG) *Y1- (Xl-ARWG) ♦ Y2) / (X2-X1) 

ENDIF 

190 CONTINUE 

CALCULATE LIFT CURVE SLOPE FOR THE CANARD. THIS ASSUMES 
A THEORETICAL 2»PI 2-D LIFT-CURVE SLOPE, AND A LINEARLY TAPERED 
CANARD PLANFORM THIS USED EQN 5.85 FROM REFERENCE 3 TO 
DETERMINE THE LIFT-CURVE SLOPE, AND EQNS 5 89 a AND b TO DETERMINE THE 
SWEEP AT THE ONE-HALF CHORD OF THE WING. 

HERE IS CALCULATED THE SWEEP AT THE ONE-HALF CHORD POINT. 

IF (SCN .NE. 0.0 .OR. SCN1 .NE. 0.0) THEN 

IF (ARCN ,EQ. 0.0) ARCN = (SPANCN1*SPANCN1) / SCN1 
IF(TRCN .EQ 0.0) TRCN = TRCN1 

TANSWP2 = TAN(SWPCN/57 . 29578) - ( (1 . -TRCN) / (1 . +TRCN) ) * (1 . / ARCN) 

A1 = 0.0 
A2 = 1.55334 
700 T 1 = TAN(Al) 

T2 = TAN(A2) 

IF(ABS (T1-TANSWP2) .LE. 0.001) THEN 
SWP2 = A1 
GO TO 750 
ENDIF 

IF(ABS (T2-TANSWP2) .LE. 0.001) THEN 
SWP2 = A2 
GO TO 750 
ENDIF 

IF (ABS (T 1 -TANSWP2) .LE. ABS(T2-TANSWP2))THEN 
IF (ABS (T 1 -TANSWP2) .ED. ABS (T2-TANSWP2)) THEN 
SWP2 = 0.785398 
GO TO 750 
ENDIF 

A2 = A2 - (A2 - All/2.0 
GO TO 700 
ELSE 

A1 = A1 * (A2 - Al)/2.0 

GO TO 700 
ENDIF 

750 CONTINUE 

HERE IS CALCULATED THE LIFT-CURVE SLOPE 


CL1 = 6.28319 * ARCN 
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c 


CL2 = (ARCN * ZMACH) **2.0 
CL3 = (ARCN / COS (SWP2) ) **2 .0 
IF (ZMACH .LE. 1.0) THEN 

CLACN = CL1 / (2.0 ♦ SQRT (CL3 
ELSE 

CLACN = CL1 / (2.0 ♦ SQRT (CL3 
ENDIF 
END IF 

CLACN = CLACN / 57.29578 


4.0 - CL2)) 
CL2 - 4.0)) 


RETURN 

END 


SUBROUTINE XACCALC (ZMACH, ZCLALFA) 

THIS SUBROUTINE IS USED TO DETERMINE THE AERODYNAMIC CENTER OF 
THE WINC THIS INCLUDES MACH NUMBER EFFECTS, WING SWEEP ANGLE, AND 
WING TAPER RATIO. THE EQUATIONS USED FOR THIS ARE TAKEN FROM 
RFFFRFNCF 3 CHAPTER 8 PAGES 484 TO 485 . THIS CALCULATES THE AC AS A 
FUNCTION OF' THE RATIO OF THE TWO-TO-THREE DIMENSIONAL L ^ F T'^ VE SL0PES 
THIS RATIO IS FOUND FROM THE CLASSIC EQUATION FOR TWO-TO-THREE 
DIMENSIONAL LIFT-CURVE SLOPES, REFERENCE 4, EQUATION 5-20. 

THIS WORKS BEST UP TO MACH OF 1 , BUT STARTS TO DEVIATE BEYOND THAT POINT 


INCLUDE ’STBCM INC ’ 

THE VARIABLE ’A' IS THE POSITION OF THE QUARTER-CHORD POINT OF 
THE WING LOCATED AT THE CENTER OF THE FUSELAGE. 


C 

C 


C 

C 

C 

C 

C 

C 


BODRAD = BDMAX / 2. 

XLECL = XLEWG - (BODRAD * TAN (SWPWG/57. 29578)) 
XTECL = XLEWG ♦ ROOTWG 
A = XLECL ♦ ((XTECL - XLECL) /4.) 


RAD = 57.29578 
PI = 3.1415926 


CLACLA = 1.0 / (1.0 - ((RAD 
XA = ((1.0 + 2 . 0*TRWG) / (1 . 0 
S TAN(SWPWC/RAD)) * CLACLA * 
XACWC = XA ♦ A 
XLE = XQCBWG - CBARWC/4 


. ZCLALFA) / (PI * ARWG) ) ) 

♦ TRWG) ) * ((SPANWG / 2.0) * 
0.3333333 


RETURN 

END 








SUBROUTINE XCGCALC(I) 
INCLUDE ’STBCM. INC’ 
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THIS SUBROUTINE IS CALLED TO CALCULATE THE RELATIVE POSITIONS 
AND MOMENTS OF EACH OF THE AIRCRAFT WEIGHTS WHICH ARE 
USED TO DETERMINE OVERALL AIRCRAFT C.G. NOTE THAT PLACEMENT 
OF MOST WEIGHTS IS ARBITRARY, BUT CAN BE ADJUSTED BY THE 
USERS WITH THE XF»** INPUTS IN /STABIN/ . 

1000 CONTINUE 

XBOD = XFBOD * BODL 
MBOD = WBODY * XBOD 
XCAN = (XqCBCN ♦ XFCAN * CBARCN) 

MCAN = WCAND * XCAN 
XHT = (XqCBHT + XFHT » CBARHT) 

MHT = WHT * XHT 
XVT = (XOCBVT ♦ XFVT * CBARVT) 

MVT = WVT * XVT 

XWING = (XQCBWG ♦ XFWG * CBARWG) 

MWG = WWING * XWING 
MAIRC = WAIRC * XLEWG 
XAPU = XFAPU * BODL 
MAPU = WAPU * XAPU 
XELT = XFELT . BODL 
MELT = WELT * XELT 
XEP = XFEP ♦ BODL 
MEP = WEP ♦ XEP 
XINST = XFINST * BODL 
MINST = WINST * XINST 
XFS = XFFS * BODL 
MFS = WFS * XFS 
XHDP = XFHDP * BODL 
MHDP = WHDP * XHDP 
XSC = XFSC * BODL 
MSC = WSC * XSC 
XPA = XFPA * BODL 
MPA = WPA * XPA 
XCARGO = XFCARGO * BODL 
MCARGO = XCARGO * WCARGO 
XCREW = XFCREW * BODL 
MCREW = WCREW * XCREW 
XBOMB = XFBOMB * BODL 
MBOMB = XBOMB * WZB 
XBB2 = XFBB2 * BODL 
MBB2 = XBB2 * WBB2 
XMISS = XFMISS * BODL 
MMISS = XMISS * WZM 
XAMMUN = XFAMMUN * BODL 
MAMMUN = XAMMUN * WZA 
XNA = XFNA * BODL 
MNA = XNA * WNA 
C 

C THE WEIGHTS MODULE DOES NOT GIVE OUT THE PIVOT WEIGHTS 

C DIRRECTLY, SO IT IS CALCULATED BY SUBTRACTING THF DIFFERENT 

C PARTS FROM THE AIRFRAME WIEGHT 

C 

WPIV = WAF-WBODY-WHT-WLG-WNA-WVT-WWJNC-WCANr> 
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t 

c 

c 


c 


1100 

c 


220 

C 

C 

C 

C 


C 

C 

c 


c 

c 

c 


c 

c 

c 


IF (WPIV .LT. 0.) WPIV = 0. 

XPIV = XFPIV * BODL 
MPIV = WPIV * XPIV 

THE SAME GOES FOR THE LIFT FAN AS GOES FOR THE PIVOTS ABOVE. 

WLFTF = WPS - WE - WFS 
IF (WLFTF .LT. 0.) WLFTF = 0. 

XLIFTF = XFLIFTF * BODL 
MLIFTF = WLFTF * XLIFTF 

WLGFRT = WFLGFRT * WLG 
XLGFRT = XFLGFRT * BODL 
MLGFRT = WLGFRT * XLGFRT 
WLGR = WLG - WLGFRT 
IF (LGFLAG .GT. 0) GO TO 1100 
XLGR = XFLGR * BODL 
CONTINUE 

MLGR = XLGR ♦ WLGR 

WZENG = WE / EN 
DO 220 J = 1 , EN 

XENG(J) = (XLEPOD(J) ♦ XFENG * PODL) 

MENG(J) = WZENG ♦ XENG(J) 

CONTINUE 

FURNISHING WEIGHT IS FOUND BY SUBTRACTING THE KNOWN WEIGHTS FROM 
THE WEIGHTS OF THE FIXED EQUIPMENT. 

WFUR = WFEQ-WAIRC-WELT-WAPU-WEP-WPA-WHDP-WINST-WSC 
IF (WFUR .LT. 0.) WFUR = 0. 

XFUR - XFFUR ♦ BODL 
MFUR = XFUR * WFUR 

PAYLOAD IS FOUND LIKE FURNISHINGS 

WPAYL = WPL-WCREW-WAMMUN-WBB2-WB0MB-WMISS-WETANK-WCARG0 
IF (WPAYL .LT. 0.) WPAYL = 0. 

XPAYL = XFPAYL * BODL 
MPAYL = XPAYL ♦ WPAYL 


CALCULATE MOMENTS FOR FUEL 
IF (FUFRAC .GT. 1.0) FUFRAOl.O 
WFWG = FUFRAC * WFP 
XWOFUEL = (XFFWG * CBARWt: - XOCPWO) 

MFWG = WFWG ♦ XWOFUEL 
WFFU3 = WFP - WFWG 
IF (WFFUS .LT. 0.0) WFFUSrrO.O 
XFUSF = (BODL * XFFFUS) 

MFFUS = WFFUS * XFUSF 

NOW INCLUDE ANY EXTERNAL TANKS DEPENDING ON THE INPUT IETANK 
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XCTIPWG = ROOTWG»TRWC 

XLETIP = XLEWG+ . 25*R00TWG+SPANWG/2 . *TAN(SWPWG/57 . 29578) 

1 - . 25*XCTIPWG 

GO TO (510,520,530,540,550,560), IETANK 
C CONDITION IETANK = 0, SINGLE TANK IN CENTERLINE OF FUS. 

XEFl = (XFEX*BODL) 

MEF = XEFl* (WZET+WZEF) 

GO TO 600 

C CONDITION IETANK = 1 , TWO TANKS ON THE SIDES OF FUS. 

510 XEFl = (XFEX*BODL) 

MEF = XEFl* (WZET+WZEF) 

GO TO 600 

C CONDITION IETANK =2, TWO TANKS AT C/4 OF WING 

520 XEFl = XQCBWC + CBARWG/4 . 

MEF = XEFl * (WZET+WZEF) 

GO TO 600 

C CONDITION IETANK =3, TWO TANKS AT WINCTIPS 

530 XEFl = (XLETIP + 0.5*XCTIPWC) 

MEF = XEFl * (WZET+WZEF) 

GO TO 600 

C CONDITION IETANK = 4, SINGLE FUS. TANK, TWO WING TANKS. 

540 WEFF = . 5*WZEF + ,5+WZET 
WEFW = WZEF + WZET - WEFF 
XEFl = (XFEX+BODL) 

XEF2 = (XQCBWG+CBARWG/4.) 

MEF = XEF 1 *WEFF+XEF2*WEFW 
GO TO 600 

C CONDITION IETANK = 5, SINGLE FUS. TANK, TWO TIP TANKS. 

550 WEFF = . 75*WZEF+ . 75+WZET 
WEFWT = WZET+WZEF-WEFF 
XEFl = (XFEX*BODL) 

XEF2 = (XLETIP+ .5+XCTIPWG) 

MEF = XEF 1*WEFF+XEF2* WEFWT 
GO TO 600 

C CONDITION IETANK = 6, TWO WING TANKS, TWO TIP TANKS. 

560 WEFW = . 75*WZEF+ . 75*WZET 
WEFWT = WZEF+ WZET -WEFW 
XEFl = (XQCBWG+CBARWG/4.) 

XEF2 = (XLETIP+ . 5*XCTIPWG) 

MEF = XEFl *WEFW+XEF2*WEFWT 
600 CONTINUE 
C 

C NOW FIND MOMENTS AND WEIGHTS OF FUEL. 

WFL = WF FUS+WFWG+ WZEF •* WZET 
MFL = MFWG+MFFUS+MFF 

DETERMINE THE SUM TOTAL OF THE WEIGHTS AND THE 
C.G. OF THE AIRCRAFT 

SUMMO = MBOD+MCAN+MHT+MNA+MVT+MWG+MAIRC+MAPU+MELT+MSC+MPTV 

1 +MEP+MINST+MHDP+MPA+MFUR+MPAYI.+MCREW+MAMM' iN-MMTSS+MCARGO 

2 +MB0MB+MFL+MLGFRT+MLGR+MLIFTF-MBB2+MREM 


C 


DO 230 J = 1 , EN 
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SUMMO = SUMMO * MENG(J) 

230 CONTINUE 


SUMWT = WB0DY+WCAND+WHT*WNA*WVT*WWING+WAIRC4-WAPU*WELT*WSC+WPIV 

1 +WEP+WINST+WHDP+WPA+WFUR*WPAYL+WZA+WZM+WLCFRT-*-WLGR*WCARGO 

2 *WCREW*WZB*WREMM-‘WFL-‘WLFTF*WE+WBB2 

XCG(I) = SUMMO / SUMrr 

THIS CALCULATES THE CENTER OF GRAVITY POSITION IN PERCENT MAC. 

CGBAR(I) = ( (XCG(I) - XLE) / CBARWG)*100. 

IF (ICGPRT . EQ . 0) GO TO 000 

THIS PRINT FLAG PRINTS THE POSITIONS AND WEIGHTS OF 

THE DIFFERENT COMPONENTS USED TO DETERMINE CORRECT PLACEMENT 

OF THE AIRCRAFTS C.G 

CALL CGPRINT 

800 CONTINUE 

THIS WRITE STATEMENT PRINTS A WARNING IF THE C.G. IS CLOSER TO THE 
MAIN GEAR THEN THE 15 DEGREE LIMIT. THIS LIMIT IS TAKEN FROM NAVY 
SPECIFICATIONS SD-24J, WHICH GIVES A GOOD DESIGN VALUE. THIS 
SPECIFICATION WAS TAKEN FROM REFERENCE 5, PAGE 40, FIGURE 4.3. 

ZMG = ZFMG . BDMAX 
TAN15 = 0.26795 
TANB = (XLGR - XCC(I)) / ZMG 
IF (TANB ,LT. TAN15) THEN 
LGFLAG = 1 

XLGR = XLGR + ( 0.01 * BODL) 

WRITE (6, 700) XCG (I) , XLGR 
GO TO 1000 
END IF 

700 FORMAT (// , 5X , 47HWARNING ! ! ! C.G. IS LESS THEN 15 DEGREES FORWARD, 

1/ , 15X, 38H0F THE MAIN GEAR IN THE VERTICAL PLANE,5X,/,9HMAIN GEAR, 
220H MOVED AFT 0 . 01 *BODL , / , 5X , 6HXCG = ,F9.5,5X,6HXMC = ,F9.5,/) 

RETURN 

END 




SUBROUTINE TOROTATE (I , IPHASE) 

THIS SUBROUTINE DETERMINES THE SIZE OF THE HORIZONTAL 
TAIL OR CANARD, NEEDED TO ROTATE THE AIRCRAFT ON TAKEOFF 
THIS USES AN ALTITUDE OF 10,000 FT FOR TAKEOFF DENSITY. 

THE METHOD USED HERE IS TAKEN FROM REFERENCE 2, CHAPTER 5 PAGE 
373 - 374. THIS IS MODIFIED TO INCLUDE A CANARD AND ANY THRUST 
OFFSET DUE TO PLACEMENT OF THE ENGINES. 
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c 

INCLUDE ’STBCM . INC ’ 

C 

C 

IF (ARHT .LE. 13.) THEN 
A = ARHT 

QCLAFHT=- . 0282* ,0359*A- ,0049*A**2* ,000345*A**3- . 000009561 *A«*4 
ELSE 

QCLAFHT = (2 . *3 . 1415926) /57 . 29578 
ENDIF 

FIRST DETERMINE THE DYNAMIC PRESSURE, AND THE LIPT "F THE WING. 

VTOR = VROT * 1077.4 
QTOR = 0.5 * .0017556 * (VT0R**2.) 

NOW CALCULATE DISTANCE PERAMETERS SPECIFYING DISTANCES FROM 
THE WING LEADING EDGE TO THE PARTICULAR COMPONENT. 

ZMG = ZFMG * BDMAX 
ZTl = ZFTl * BDMAX 
ZT2 = ZFT2 * BDMAX 
ZD = ZFD * BDMAX 
T1RAD = GAMAT1 / 57.29578 
T2RAD = GAMAT2 / 57.29578 
XLEMG = XLGR - XQCBWG 
XLECG = XCG(I) - XQCBWG 
XLETl = XqCBWG - (XT1*B0DL) 

XLET2 = (XT2.B0DL) - XQCBWG 
XLEACWG = XACWG - XQCBWG 
IF (XLECN .EQ. 0.0) THEN 
XLECN = XQCBWG - XQCBCN 
ENDIF 

XLEHT = XqCBHT - XQCBWG 

CALCULATE THE LIFT OF THE WING. THIS INCLUDES GROUND EFFECTS 
FOR THE LIFT COEFFICIENT. THIS IS JUST MULTIPLICATION OF THE 
WING LIFT-CURVE SLOPE MODIFIED FOR GROUND EFFECT, 

BY THE ANGLE OF ATTACK OF THE WING, AND THE DYNAMIC PRESSURE 
FIRST MODIFY FOR GROUND EFFECTS 

DG = (ZRTWG+ZMG) / (SPANWC/2.) 

CALL GRNDEFFECT (DG , AFGDAF) 

CLAGRD = AFGDAF * CLALFA 
LWB = CLAGRD * (ALFTO) * QTRR * "WO 

CALCULATE THE CHANGE IN LIFT-CURVE SLOPE OF THE TAIL WHILE 
IN GROUND EFFECT. 

IF (SPANHT . LE. 0.0) GO TO 13 
DG = (ZRTHT+ZMG) / (SPANHT/2.) 

CALL GRNDEFFECT (DG, AFGDAF) 

CLAHTGRD = AFGDAF * QCLAFHT 
13 CONTINUE 
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c 

c 

C CALCULATE THE CHANCE IN LIFT-CURVE SLOPE OF THE CANARD WHILE 
C IN GROUND EFFECT. 

C 

IF (SPANCN .LE. 0.0) GO TO 14 
DG = (ZRTCN+ZMG) / (SPANCN/2.) 

CALL GRNDEFFECT (DG , AFGDAF) 

CLACNGRD = AFGDAF . CLACN 
14 CONTINUE 

IF (SPANCN1 .LE. 0.0) GO TO 41 
DG = (ZRTCN1+ZMG) / (SPANCNl/2.) 

CALL GRNDEFFECT (DC, AFGDAF) 

CLACNGRD = AFGDAF * CLACN 
41 CONTINUE 
C 

C NOW CALCULATE THE DOWNWASH ANGLE OF THE WING ON THE TAIL. 

C 

IF (SHT .EG. 0.0) GO TO 10 
CALL DOWNWASH (CLTO,ALFTO,DWANGLE) 

10 CONTINUE 
C 

C THE CALL TO TRAJDT IS USED TO CALCULATE THE MAXIMUM THRUST AT 

C TAKEOFF ROTATION OF THE AIRCRAFT. 

C 

M0DN04=4 

CALL TRAJDT (M0DN04 , ICALC , NERROR , IGEO , KGPRNT , IGPLT) 
IF(NERR0R.GE.2) RETURN 
TN = THRUST *EN 

C 

TWZT0=TN 

C 

C NOW CALCULATE THE THRUST SPLIT, IF ANY. 

SHTTEMP = SHTl 
SCNTEMP = SCN1 
SHTl = SHT 
SCN1 = SCN 

IF (IVECT .NE. 0 .AND. TSPLIT .EQ. 0.0) THEN 
CALL VECTHRUST (I) 

END IF 

SHTl = SHTTEMP 
SCN1 = SCNTEMP 
TJ = TSPLIT + TWZTO 
T2 = TWZTO - Tl 

NOW CALCULATE THE LIFT OF THE HORIZONTAL TAIL AND CANARD 
DEPENDING ON WHICH OF THE TWO EXIST. IF BOTH EXIST, THE 
CANARD IS ASSUMED TO BE USED ONLY FOR TRIM, THE HORIZONTAL TAIL 
IS THE DESIRED SIZE. 

THE METHOD USED HERE IS TAKEN FROM REFERENCE 2, CHAPTER 5 PAGE 
373 - 374. THIS IS MODIFIED TO INCLUDE A CANARD AND ANY THRUST 
OFFSET DUE TO PLACEMENT OF THE ENGINES. 

IF (SCN1 .EC). 0.0 .AND. SCN .EQ 0.0) THEN 
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C ** THIS IS THE CASE WHERE THERE IS A HORIZONTAL TAIL ONLY * 
LH1 = CMWBTO * QTOR * SWG * CBARWG 
LH3 = ZTl * T1 
LH5 = ZT2 * T2 

LH6 = WGTO * (XLEMG - XLECG + ZMU * ZMG) 

LH7 = LWB * (XLEMG - XLEACWG * ZMU * ZMG) 

LH8 = D * ZD 

DENOM = XLEHT - XLEMG - ZMU * ZMG 

LHT = (LH1 + LH2 + LH3 - LH4 * LH5 - LH6 * LH7 - LH8) / 

C 

SHTl=LHT / (CLAHTGRD*QTOR* (ALFTO-DWANGLE*IHT*TAU*ELVDEF)) 
ELSE 

IF (SHT .EQ. 0.) THEN 

•* THIS IS THE CASE WHERE THERE IS A CANARD BUT NO H. TAIL 

LC1 = CMWBTO * QTOR » SWG * CBARWG 
LC3 = ZTl * Tl 

LC5 = ZT2 * T2 

LC6 = WGTO * (XLEMG - XLECG * ZMU * ZMC) 

LC7 = LWB * (XLEACWG - XLEMG - ZMU * ZMG) 

LC8 = D * ZD 

DENOM = XLECN * XLEMG - ZMU * ZMG 

LCN = (-LC1 - LC2 - LC3 * LC4 - LC5 ♦ LC6 ♦ LC7 ♦ LC8) / 

SCN1 = LCN/ (CLACNCRD*QTOR* (ALFTO+ICN+TAU*ELVDEF) ) 

ELSE 

*• THIS IS THE CASE WITH BOTH CANARD AND HORIZONTAL TAIL 

LH1 = CMWBTO * QTOR * SWG * CBARWG 
LH3 = ZTl • Tl 

LH5 = ZT2 * T2 

LH6 = WGTO * (XLEMG - XLECG + ZMU * ZMG) 

LH7 = LWB * (XLEMG - XLEACWG «■ ZMU ♦ ZMG) 

LH8 = D * ZD 

LCN = CLACNGRD * QTOR * SCN1 * (ALFTO ♦ ICN) 

LH9 = LCN * (XLECN ♦ XLEMG + ZMU * ZMG) 

DENOM = XLEHT - XLEMG - ZMU * ZMG 

LHT = (LH1 + LH2 * LH3 - LH4 ♦ LH5 - LH6 - LH7 - LH8 
J ♦ LH9) / DENOM 

C 

SHT1=LHT/ (CLAHTGRD'QTOR. (ALFTO-DWANCLE-JHT-TAU'ELVDFF') ) 
ENDIF 
ENDIF 

C 

C 

RETURN 

END 

C 

f ♦♦•*******#*******************i»’****#*‘********'*'*’ + ’*-*'*i**»*:*******' 

C 

SUBROUTINE XCGFORWARD (I) 

C 


* 


DENOM 


* * 


DENOM 


* * 
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THIS SUBROUTINE CALCULATES THE LONCITUDINAL CONTROL SURFACE 
SIZE NEEDED TO MAINTAIN CONTROLABILITY OF THE AIRCRAFT IN 
LANDING CONFIGURATION. THAT IS MAX CL, LOW SPEED, HIGH DRAG, 

IN GROUND EFFECT. THIS IS DONE BY A SUMMATION OF MOMENTS ABOUT 
THE CENTER OF GRAVITY. WITH THIS EQUATION SET TO ZERO, THE 
CENTER OF GRAVITY POSITION IS SOLVED FOR. 


INCLUDE ’STBCM . INC ’ 

VTHT = (CLAHTGRD * SHT1 * ETAHT)/ SWG 
VTCN = (CLACNGRD * SCN1 * ETACN) / SWG 

IF (SHT .EQ. 0.0) GO TO 20 

CALL DOWNWASH (CLMAX , ALFMAX , DWANGLE) 

20 CONTINUE 

DETERMINE THE DIFFERENT ANGLES OF ATTACK. 

ALFHT = ALFMAX •* IHT - DWANGLE 
ALFCN = ALFMAX ♦ ICN 
IF (SCN .EQ. 0.0) THEN 

ALFHT = ALFHT - TAU * ELVDMAX 
ELSE 

IF (SHT .EQ. 0.0) THEN 

ALFCN = ALFCN + TAU * ELVDMAX 
ALFHT =0.0 
ELSE 

ALFHT = ALFHT + TAU * ELVDMAX 
END IF 
ENDIF 


XAWBZ = XLE + XAWB * CBARWG 
CM1 = CLMAX * XAWBZ 
CM2 = CMWBMAX * CBARWG 
CM3 = VTHT * ALFHT . XQCBHT 
CM4 = VTCN * ALFCN « XQCBCN 
CM5B = ZT1 * T1 
CM5D = ZT2 * T2 
CM5 = (CM5B ♦ CM5D) * CBARWG 

DENOCM = CLMAX ♦ (VTHT * ALFHT) * (VTCN » ALFCN) 
THIS SOLVES FOR THE FORWARD C.G. LIMIT. 


XCGFRWD (I) = (CM1 - CM2 * CM? - CM4 - CM5) / nFNOCM 
C 

IF (IDBPRT . NE . 0)THEN 

WRITE (6 , 10) 1-2 , XCGFRWD (T ) 

ENDIF 

C 

10 FORMAT (10X,8HXCGFRWD[, 12, 4H) = ,F12.6) 

RETURN 

END 

C 

£ ++*+**+++***++++++++++*++**+***++**+++*+*+*+***+********** 


******** 
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SUBROUTINE VECTHRUST(I) 

THIS SUBROUTINE SOLVES FOR ONE OF THREE SPECIFIED VARIABLES IN THE 
SUMMATION OF FORCES ABOUT THE C.G., THE THRUST ANGLE OF THE REAR 
NOZZLE THE POSTION OF THE FOREWARD THRUST NOZZLE, AND THE THRUST 
SPLIT OF THE FORE AND AFT NOZZLES. THE VALUES CALCULATED ARE THOSE 
NEEDED TO EQUATE THE MOMENTS ABOUT THE C.G. TO ZERO. 


INCLUDE ’STBCM . INC ’ 

WEMP=WCT0-WFT0T-WMISS-WAMMUN-WB0MB-WBB2-WBB1 -WETANK 
Vl = CLMAX . ( (XCG(I) - XACWG) / CBARWG) 

V2 = CMWBMAX 

V3A = CLAFHT * ALFHT * ETAHT * SHT1 / SWG 
V3B = (XQCBHT - XCC(I)) / CBARWG 
V3 = V3A * V3B 

V4A = CLACN * ALFCN . ETACN * SCN1 / SWG 
V4B = (XCG(I) - XQCBCN) / CBARWG 
V4 = V4A * V4B 
SUMVT = VI + V2 - V3 + V4 


C 


100 


110 


120 


150 


GO TO (100,110,120), IVECT 
CONTINUE 

IF(T1*SIN(T1RAD) ♦ T2*SIN (T2RAD) .LT. WEMP) THEN 
WRITE (6 , 500) T1 *SIN (T1RAD) +T2»SIN(T2RAD) , WEMP 
GO TO 170 
ENDIF 

V5 = ZT1 * Tl * COS (T1RAD) 

V6 = ((XT 2 * BODL) - XCG(I)) * T2 * SIN(T2RAD) 

V7 = ZT2 * T2 * COS (T2RAD) 

SUMVTT = SUMVT ♦ V5 - V6 + V7 

XTl = (SUMVTT / (Tl * SIN(TIRAD)) * XCG(I)) / BODL 

GO TO 170 

CONTINUE 

NUMVT = (SUMVT /TWZTO) - ( (XT2*B0DL) -XCG (I) ) *SIN (T2RAD) 

$ +ZT2*C0S(T2RAD) 

DENOVT = ZT2*C0S (T2RAD) - ( (XT2*B0DL) -XCG (I) ) *SIN (T2RAD) 

J - (XCG (I) - (XT1*B0DL) ) *SIN (T1RAD) -ZT1*C0S (TlRAD) 

TSPLIT = NUMVT / DENOVT 

GO TO 170 

CONTINUE 

DUMM=SUMVT+ (XCG(I) - (XT1*B0DL) ) *T 1*SIN(T1RAD) *7Ti -ti *COS (T1RAP) 
GAM21 =0.0 
GAM22 = 1.5707963 


DUM1 = DUMM+ZT2»T2*C0S (GAM21) - (XT2*B0DL-XCG(I)) - T 2*SIN(GAM21) 

DUM2 = DUMM+ZT2*T2*C0S(CAM22)-(XT2*B0DL-XCG(T)W2*SIN(GAM22) 

IF (ABS (ABS (DUM1 ) -ABS (DUM2) ) LT. 0.05)THEN 
T2RAD = GAM21 
GO TO 170 
ENDIF 

IF (ABS(DUMl) .LT. ABS(DUM2)) THEN 

GAM22 = GAM22 - (GAM22 - GAM21) / 2.0 





ELSE 

GAM21 = GAM21 * (GAM22 - GAM21) / 2.0 
ENDIF 
GO TO 150 
170 CONTINUE 
C 

500 FORMAT (10X, 61HWARNIG ! ! VERTICAL THRUST COMPONENT LESS THEN 
J AIRCRAFT WEIGHT ! , /, 10X, 7HTHRUST=, F15 .7,5X,7HWEIGHT=, F15 . 7) 
RETURN 
END 




SUBROUTINE DATATRANS 
INCLUDE 'STBCM . INC ' 

THIS SUBROUTINE IS USED TO CALL THE DATA TRANSFER SUBROUTINE. IT 
IS USED TO DETERMINE THE ANGLE OF ATTACK, COEFFICIENT OF LIFT, AND 
PITCHING MOMENT COEFFICIENT FOR TAKEOFF, MAX LIFT CONDITIONS, 

AND ZERO ANGLE OF ATTACK CONDITIONS. 

MACH =0.35 
ALT = 5000.0 
CL = 0.54 
CD = .54 
CMWB =0.54 
ALPHA = AWTOT 
ICALC = 2 
IAO = 1 

CALL STBDT (ICALC , NERR , IGEO , KGPRNT , IGPLT) 

ALFTO = ALPHA 
CLTO = CL 
CDTO = CD 
CMWBTO = CMWB 
C 77 . 7 , 7 . 7 . 

CL = 0.54 
CD = .54 
CMWB =0.54 
ALPHA =0.54 
IAO = 10 

CALL STBDT (ICALC , NERR , IGEO , KGPRNT , TGPLT) 

ALFMAX = ALPHA 
CLMAX = CL 
CDMAX = CD 
CMWBMAX = CMWB 
C 77 . 7 . 7 , 7 . 

CL = 0.0001 
CD = .54 
CMWB =0.54 
ALPHA =0.54 
IAO = 0 

CALL STBDT (ICALC , NERR , IGEO , KGPRNT , IGPLT) 

ALFO = ALPHA 
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CLO = CL 
CDO = CO 
CMWBO = CMWB 

RETURN 

END 




SUBROUTINE DOWNWASH (CLZ , ALFZ, DWANGLE) 

THIS SUBROUTINE CALCULATES THE DOWNWASH ANGLE ON THE TAIL. FIRST 
BY THE THEORETICAL METHOD, AND THEN BY POSITION OF THE TAIL RELATIVE 
TO THE WING VOTECIES. 

INCLUDE ’STBCM . INC ’ 

CALCULATE THE DOWNWASH ANGLE OF THE WING ON THE TAIL. 

DWANGLE = ((1.62 * CLZ) / (3.1415926 * ARWG) ) * 57.29578 

NOW MODIFY THIS CALCULATION FOR THE POSITION OF THE TAIL. THE 
METHOD USED IS FROM REFERENCE 1, SECTION 4.4.1. 

FIRST IS A CURVE FIT OF FIGURE 4.4.1-55. 

ALFRATIO = (ALFZ - ALFO) / (ALFMAX - ALFO) 

IF (SWPWG LE. O.O) THEN 
AREFF = ARWG 
SPANEFF = SPANWG 
GOTO 1100 
ENDIF 

IF (ALFRATIO .LT. 0.4) THEN 
AREFF = ARWG 
SPANEFF = SPANWG 
GOTO 1100 
ENDIF 

IF (ALFRATIO .LT. 0.56 .AND. SWPWG .LT. 45.0) THEN 
AREFF = ARWG 
SPANEFF = SPANWG 
GOTO 1100 
ENDIF 

IF (ALFRATIO .LT. 0.56) THEN 

A1 = -2.5098 «■ 6.2131 «ALFRATIP 
GOTO 1000 
ENDJF 

IF (ALFRATIO .LT. 0.73 AND. SWFWO LT. 30.0) THFM 
AREFF = ARWG 
SPANEFF = SPANWG 
GOTO 1100 
ENDIF 

IF (ALFRATIO .LT. 0.73) THEN 
XI = 45.0 
X2 = 60.0 

Y1 = -2.8585 ♦ 5 . 0852+ALFRATIO 
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Y2 = -2.5098 + 6 . 2131*ALFRATI0 

Al = ( (X2-SWPWG) *Y1- (Xl-SWPWG) +Y2) / (X2-X1) 

IF (A1 . LT . 0.0) Al = 0.0 
GOTO 1000 
ENDIF 

IF (ALFRATIO .CE. 0.73) THEN 
IF (SWPWG . LE. 45.0) THEN 
XI = 30.0 

X2 = 45.0 

Y1 = -2.6464 + 3 . 6316*ALFRATI0 

Y2 = -2.5098 ♦ 6 . 2131 *ALFRATIO 

Al = ((X2-SWPWG) *Y1- (Xl-SWPWG) *Y2)/ (X2-X1) 

IF (Al .LT. 0.0) Al = 0.0 
ELSE 

XI = 45.0 

X2 = 60.0 

Y1 = -2.8585 + 5 ,O852*ALFRATI0 

Y2 = -2.5098 + 6.2131*ALFRATI0 

Al = ((X2-SWPWG)*Y1- (Xl-SWPWG) *Y2)/(X2-X1) 

IF (Al .LT. 0.0) Al = 0.0 
ENDIF 
ENDIF 

1000 CONTINUE 

IF (TRWG .LT. 0.50) THEN 
XI = 0.0 
X2 = 0.5 

YJ = 1 .0 - 0 . 150*A1 
Y2 = 1.0 - . 1245*A1 

AEFA = ((X2 - TRWG) *Y1- (Xl-TRWG) *Y2) / (X2-X1) 

AREFF = AEFA * ARWG 
ELSE 

XI = 0.5 
X2 = 1.0 

Y1 = 1.0 - . 1245*A1 
Y2 = 1.0 - . 1000*A1 

AEFA = ( (X2-TRWG) *Y1 - (XI -TRWG) »Y2) / (X2-X1) 

AREFF = AEFA * ARWG 
ENDIF 
XI = 0.0 
X2 = 10 

Y1 = - ! 20326+2 . 6606*AEFA-2 . 2631* (AEFA* *2 .)-*■■ 80556* (AEFA* *3 ) 
Y2 = -.016464 + 1.0168*AEFA 
SPANEFB = ( (X2-TRWG) *Y1- (Xl-TRWG) »Y2) / (X2-X1I 
SPANEFF = SPANEFB * SPANWO 
1100 CONTINUE 

NOW SOLVE FOR THE VERTICAL POSITION OF THE VORTEX CORE. 

THIS ASSUMES TRAILING EDGE SEPERATION . THESE USE EQUATIONS 
4.4.1-A THRUOGH 4.4.1-F IN REFERENCE 1, SECTION 4.4.1. 


SICRU = (0.56 ♦ ARWG) / CL7 

BVRU = (0.78 + 0.10*(TRWG-0.4)*0.003*SWPWG)*SPANEFF 
LEFF = XQCBHT - XLETIP 

BV = SPANFFF - (SPANEFF - BVRU) -SORT ( (2 O-LEFF) / (SPANWC- 


?TGRU)) 



nonnnn r*i o o o r~) ri 


99 


HH = ZRTHT - ZRTWG 

HVERT = HH - LEFF» ( (ALFZ/57 . 29578) - (0.41*CLZ)/(3.14159*AREFF)) 
$ - (SPANEFF/2 . 0) *TAN (OTHED/57 . 29578) 

NOW CALCULATE THE MODIFIED DOWNWASH ANGLE. 

DWANGLE = DWANGLE * (1 / (1 * ((2,0*HVERT)/BV)..2,)) 

RETURN 

END 




SUBROUTINE GRNDEFFECT (DG , AFGDAF) 

THIS SUBROUTINE DETERMINES THE CHANGES TO CLALPHA FOR TAKEOFF 
AND LANDING CALCULATIONS. THIS USES FIGURE 5-39 FROM REFERENCE 
4, PG. 257. THE HEIGHT FROM THE GROUND USED IS THE HEIGHT OF 
AIRCRAFT WITH THE MAIN LANDING GEAR ON THE GROUND. 

INCLUDE ’STBCM . INC ’ 

C 

IF (ARWG .LE. 6.) THEN 
XI = 4. 

X2 = 6 

Yl=l .2989- . 72851 +DG* . 79231 *DC* *2 . - . 4131*DG»»3 . + .083452*DG»*4 . 
Y2=l .2054- .49627.DG* . 54742*DG**2 . - . 29063*DG**3 . - ,059777*DG**4 . 
AFCDAF=( (X2-ARWG) *Y1- (Xl-ARWG) .Y2) /(X2-X1) 

GO TO 190 
END IF 

IF (ARWG ,GT. 6. AND. ARWG .LE. 8.) THEN 
XI = 6. 

X2 = 8. 

Yl=l . 2054- . 49627.DG+ . 54742*DG**2.-. 29063*DG**3 . + .059777.DG**4 
Y2=l . 1751- . 46048.DG+ . 55655.DG..2 . - . 31668*DC**3 . + .068063.DC+.4 
AFGDAF= ( (X2-ARWG) * Y1 - (XI -ARWG) * Y2) / (X2-X1) 

GO TO 190 
ENDIF 

IF (ARWG . GT . 8. .AND. ARWG .LE. 10.) THEN 
XI = 8. 

X2 = 10. 

Yl=l . 1751- . 46048*DC+ . 55655*DG**2 . - . 31668*DG*-3 .068063*00*. 4 
Y2=l . 1498- 41037.DG* . 51?26*DG*»2 . - . 30114.DG**? . * ,06628«*l»i;**/i 
AFGDAF= ( (X2-ARWG) » Yl- (Xl -ARWG) *Y2) / (X2-X1 ) 

GO TO 190 
ENDIF 

IF (ARWG . GT . 10.) THEN 
XI = 10. 

X2 = 12. 

Y1=J .1498- .41037.DG* . 51326»DG*»2 . - . 30114.DG* *3 . * . 066288.DG* .4 
Y2=l .1188- . 31846* DG* .3826* DGt*2. - .21 725*DG*»? * .046757*06* *4 
AFGDAF= ( (X2-ARWG) .Yl - (XI -ARWG) * Y2) / (X2-X1 ) 

ENDIF 
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190 CONTINUE 

c 

RETURN 

END 




SUBROUTINE STABOUT 

THIS SUBROUTINE PRINTS OUT THE PERTINATE OATA CALCULATED BY THE 
MODULE . 

INCLUDE ’STBCM . INC ’ 

C 

DO 53 I = 3.NPHASE-3 

IF (XCGFRWD (I) .CT. XCGFRWD(1))THEN 
XCGFRWD (1) = XCGFRWD (I) 

END IF 

53 CONTINUE 
C 

DO 54 I = 3.NPHASE+3 

IF (XCGAFT (I) LT. XCGAFT (2) )THEN 
XCGAFT (2) = XCGAFT (I) 

END IF 

54 CONTINUE 
C 

DO 55 I = 2 ,NPHASE+2 

SM(I) = ( (XCGAFT (2) - XCG(I))/CBARWG).100. 

55 CONTINUE 
C 

WRITE (6, 50) 

WRITE (6, 105) (XAWB*CBARWG) +XLE 
IF (SHT1 .NE. 0.) THEN 
WRITE(6, 150) SHT1 
ENDIF 

IF (SCN1 .NE. 0.) THEN 
WRITE (6, 160) SCN1 
ENDIF 

WRITE(6, 135) XCGFRWD (1) 

WRITE (6, 140) XCGAFT (2) 

DO 10 I = 2,NPHASE*2 
IF (I . LE. 2) THEN 
WRITE(6,59) 

ELSE 

WRITE (6, 60) I - 2 
ENDIF 

WRITE (6, 120) FNEUTPT (I) 

WRITE(6, 100) XCG(I) 

WRITE(6,107) CGBAR(I) 

WRITE (6, 110) DCMDCL (I) 

WRITE (6, 130) SM(I) 

10 CONTINUE 
C 

IF (IVECT ,EQ. 0) THEN 
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GO TO 500 
ELSE 

WRITE (6,1 70)T1*SIN(T1 RAD) , T2*SIN(T2RAD) 
WRITE (6, 180) XT 1 ,XT2 
WRITE(6, 185) ZTl , ZT2 

WRITE (6 , 190) T1 RAD* 57 . 29578 ,T2RAD*57 . 29578 
WRITE(6,195)TSPLIT 
END IF 

500 CONTINUE 

C 


WRITE (6, 200) 


50 

59 

60 
100 
105 
107 
110 
120 
130 
135 
140 
150 
160 
170 

180 

185 

190 

195 

200 


FORMAT (1H1 , 10X , 21HSTABILITY OUTPUT DATA,//) 
FORMAT (5X,23HDATA AT END OF TAKE-OFF) 

FORMAT (5X,2OHMISSI0N PHASE NUMBER, 15) 


FORMAT (10X,12HXCG 
FORMAT (10X,12HXAC 
FORMAT (10X,12HCGBAR = 
FORMAT (10X,12HDCMDCL = 
FORMAT (10X.12HNEUTPT 


, F15 . 7 ,4H ft.) 

, F15 . 7,4H ft.) 

, F15.7,8H % chord) 

, F15 . 7) 

, F15 . 7) 

FORMAT (10X, 16HSTATIC MARGIN = ,F15.7,8H % chord) 

FORMAT (1 OX, 12HXCGFRWD = ,F15.7,4H ft.) 

FORMAT (10X, 12HAFT CG = ,F15.7,4H ft.) 

FORMAT (10X, 12HTAIL SIZE = ,F15.7,8H sq. ft.) 

FORMAT (10X, 14HCANARD SIZE = ,F15.7,8H sq. ft.) 

FORMAT (10X,15HVERTICAL THRUST, /,10X,12HFORE THRUST=,F15.7, 
S5X, 12HAFT THRUST =,F15.7) 

FORMAT (10X,21HFORE THRUST P0SITI0N=,F15.7, 

J5X,20HAFT THRUST P0SITI0N=,F15.7) 

FORMAT (1 OX, 22HVERT FORE THRUST P0S.=,F15.7,5X, 

S22HVERT AFT THRUST POS. =,F15.7) 

FORMAT (1OX.21HF0RE THRUST ANGLE = ,F15.7,5X, 

J22HAFT THRUST ANGLE = ,F15.7) 

FORMAT (10X,24HTHRUST SPLIT (FORE/TOT) =,F15. 7) 

FORMAT (/,10X,20HEND STABILITY OUTPUT,//) 


RETURN 

END 




SUBROUTINE CGPRINT 
INCLUDE ’STBCM. INC ' 

THIS SUBROUTINE PRINTS OUT THE COMPONENT WEIGHTS AND 
POSITIONS WHEN ICGPRT = 1 . THIS IS USED TO DETERMINE 
THE CORRECT POSITION OF THE C.C. OF THE AIRCRAFT. 

WRITE (6,700) 

WRITE (6,702) 

WRITE (6,705) W80DY,XB0D 
WRITE (6,710) WWING , XWINC 
WRITE (6,715) WCAND, XCAN 
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WRITE (6,720) WHT , XHT 
WRITE (6,725) WVT,XVT 
WRITE (6,730) WLCFRT,XLCFRT 
WRITE (6,732) WLGR ,XLGR 
WRITE (6,735) WNA , XNA 
WRITE (6,740) WPIV,XPIV 
WRITE (6,745) WAIRC,XLEWC 
WRITE (6,750) WAPU,XAPU 
WRITE (6,760) WELT,XELT 
WRITE (6,765) WEP,XEP 
WRITE (6,770) WIN5T,XINST 
WRITE (6,775) WHDP,XHDP 
WRITE (6,780) WSC,XSC 
WRITE (6,785) WPA,XPA 
WRITE (6,787) WFUR,XFUR 
DO 795 J = 1,EN 
WRITE (6,805) WZENG,XENG(J) 

795 CONTINUE 

WRITE (6,810) WFS,XFS 
WRITE (6,815) WLFTF , XLIFTF 
WRITE (6,820) WPAYL,XPAYL 
WRITE (6,830) WCREW,XCREW 
WRITE (6,835) WCARGO ,XCARCO 
WRITE (6,840) WAMMUN , XAMMUN 
WRITE (6,855) WBB2,XBB2 
WRITE (6,860) WBOMB , XBOMB 
WRITE (6,865) WMISS,XMISS 
WRITE (6,880) WFFUS,XFUSF 
WRITE (6,890) WFWG , XWGFUEL 
WRITE (6,870) 

IF (IETANK .EQ. 0) THEN 
WRITE (6,1050) 

WRITE (6,1051) XEF1 
WRITE (6,1052) WZET,WZEF 
ELSE 

GOTO (871,872,873,874,875,876), IETANK 

871 WRITE (6,1100) 

WRITE (6,1101) XEF1 
WRITE (6,1102) WZET, WZEF 
GO TO 900 

872 WRITE (6,1200) 

WRITE (6,1201) XEFl 
WRITE (6,1202) WZET, WZEF 
GO TO 900 

873 WRITE (6,1300) 

WRITE (6,1301) XEFl 
WRITE (6,1302) WZET, WZEF 
GO TO 900 

874 WRITE (6,1400) 

WRITE (6,1401) XEFl , XEF2 
WRITE (6,1402) WEFF.WEFW 
GO TO 900 

875 WRITE (6,1500) 

WRITE (6,1501) XEFl ,XEF2 
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WRITE (6,1502) WEFF , WEFWT 
GO TO 900 

876 WRITE (6,1600) 

WRITE (6,1601) XEF1.XEF2 
WRITE (6,1602) WEFW, WEFWT 
GO TO 900 

900 CONTINUE 
END IF 

WRITE (6,1700) 

C 

700 FORMAT (/,5X,40HAIRCRAFT COMPONENT WEICHTS AND POSITIONS,/) 
702 FORMAT (5X, 12H COMPONENT ,12H WEIGHT ,12H POSITION ) 
705 FORMAT (5X, 12H AIRFRAME ,F12 ,6,F12.6) 

710 FORMAT (5X, 12H WING ,F12.6,F12.6) 

715 FORMAT (5X , 12H CANARD ,F12.6,F12.6) 

720 FORMAT (5X, 12H HT ,F12 6,F12.6) 

725 FORMAT (5X , 12H VT , F12 . 6 , F12 .6) 

730 FORMAT (5X, 12H NOSE GEAR ,F12.6,F12.6) 

732 FORMAT (5X, 12H MAIN GEAR ,F12.6,F12.6) 

735 FORMAT (5X,12H NACELLS ,F12.6,F12.6) 

740 FORMAT (5X, 12H PIVOTS ,F12.6,F12.6) 

745 FORMAT (5X, 12H AIR COND. ,F12.6,F12.6) 

750 FORMAT (5X, 12H APU , F12 . 6 , F12 . 6) 

760 FORMAT (5X, 12H AVIONICS , F12 . 6 , F12 . 6) 

765 FORMAT (5X, 12H ELECTRICAL, F12. 6, F12. 6) 

770 FORMAT (5X,12H INSTURMENTS,F12.6,F12.6) 

775 FORMAT (5X, 12H HYDRAULICS ,F12.6,F12.6) 

780 FORMAT (5X, 12HC0NTR0L SURF,F12.6,F12.6) 

785 FORMAT (5X, 12H PA ,F12.6,F12.6) 

787 FORMAT (5X, 12H FURNISHINGS, F12. 6, F12. 6) 

805 FORMAT (5X,12H ENGINES ,F12.6,F12.6) 

810 FORMAT (5X, 12H FUEL SYS. ,F12.6,F12.6) 

815 FORMAT (5X, 12H LIFT FAN ,F12.6,F12.6) 

820 FORMAT (5X, 12H PAYLOAD ,F12.6,F12.6) 

830 FORMAT (5X, 12H CREW , F12 . 6, F12 .6) 

835 FORMAT (5X, 12H CARGO ,F12.6,F12.6) 

840 FORMAT (5X,12H AMMUNITION ,F12.6,F12.6) 

855 FORMAT (5X, 12H BB2 , F12 . 6, F12 .6) 

860 FORMAT (5X, 12H BOMBS , F12 . 6, F12 .6) 

865 FORMAT (5X,12H MISSILES ,F12.6,F12.6) 

870 FORMAT (5X, 25HEXTERNAL TANK (S) AND FUEL) 

880 FORMAT (5X,12H FUS . FUEL ,F12.3,F12.6) 

890 FORMAT (5X,12H WING FUEL ,F).2.3,F12.6) 

1050 FORMAT (5X , 20HCENTERLINE FUEL TANK) 

1051 FORMAT (5X,15HEXT. TANK P03.=,F12.6) 

1052 FORMAT (5X,8HW. TANK= , F 1 2 . 6 , 5X , 8HW FUEL=,F12.6) 

1100 FORMAT (5X,19HTW0 FUS. SIDE TANKS) 

1101 FORMAT (5X , 15HEXT . TANK P0S.=,F12.6) 

1102 FORMAT (5X,8HW. T ANK= , F 1 2 . 6 , 5X , 8HW FUEL=,F12.6) 

1200 FORMAT (5X,18HTW0 C/4-WTNG TANKS) 

1201 FORMAT (5X,15HEXT. TANK P0S.=,F12.6) 

1202 FORMAT (5X,8HW. T ANK= , F 1 2 . 6 , 5X , 8HW FUEL=,F12 *) 

1300 FORMAT (5X,17HTW0 WINGTIP TANKS) 
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1301 FORMAT (5X, 

1302 FORMAT (5X r 

1400 FORMAT (5X, 

1401 FORMAT (5X, 

1402 FORMAT (5X, 

1500 FORMAT (5X, 

1501 FORMAT (5X, 

1502 FORMAT (5X, 

1600 FORMAT (5X, 

1601 FORMAT (5X, 
S18HWINGTIP 

1602 FORMAT (5X, 
1700 FORMAT (//) 

RETURN 

END 


15HEXT . TANK P0S.=,F12.6) 

8HW . TANK=, F12 . 6 , 5X , 8HW . FUEL=,F12.6) 

36H0NE C.L. TANK AND TWO C/4-WING TANKS) _ 

14HC.L TANK P0S.=,F12.6,5X,19HC/4-WING TANK P0S.-,F12.6) 
13HW FUS FUEL=,F12.6,5X,13HW. WING FUEL=,F12.6) 

35H0NE C.L. TANK AND TWO WINCTIP TANKS) 

15HC.L. TANK POS .=, F12 . 6 ,5X, 18HWINGTIP TANK P0S.=,F12.6) 
13HW. FUS. FUEL=,F12 .6 ,5X,13HW. WING FUEL=,F12.6) 

37HTW0 C /4-WING TANKS, TWO WINGTIP TANKS) 

19HC/4-WING TANK P0S.=,F12.6,5X, 

TANK POS .=, F12 . 6) 

13HW . WING FUEl=,F12.6,5X,16HW. WINGTIP FUEL=,F12.6) 






************************* 


SUBROUTINE STBDT (ICALC , NERR , ICEO , KGPRNT , IGPLT) 

THIS SUBROUTINE IS USED TO TRANSFER DATA BETWEEN THE STABILITY 
MODULE AND THE AERODYNAMICS MODULE. IT IS IMPORTANT TO MATCH 
THE GLOBAL ARRAY SIZES BETWEEN THIS SUBROUTINE AND THE MODULE 
’DATATR.DAT’. NOTE THAT FOR EACH COMMON BLOCK OF THE MODULES 
THERE EXISTS ONE EXTRA VARIABLE 

THEN AS CALLED OUT IN DATATR.DAT. THIS ALLOWS FOR A DUMMY 
VARIABLE TO BE INCLUDED. 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 


C 

C 


COMMON/GLOBCM/RARAY (1400) , IARAY (100) 
C0MM0N/ACSNT2/IN0UT (1600) , ILOCAT (30,4) 
STABILITY AND CONTROL COMMON -MODULE 5. 

COMMON /STBCM/ RA5(134) ,IA5(2) 
AERODYNAMICS COMMON - MODULE 3. 

COMMON /AEROCM/ RA3(170) ,IA3(10) 


DATA TRANSFER INTERFACE 
OPERATION SEQUENCE. 

1 - TRANSFER DATA FROM 

2 - TRANSFER DATA FROM 

3 - CALL AERODYNAMICS. 

4 - TRANSFER DATA FROM 

5 - TRANSFER DATA FROM 


BETWEEN STABILITY AND AERODYNAMICS. 

STABILITY LOCAL COMMON TO GLOBAL COMMON. 
GLOBAL TO AERODYNAMICS LOCAL COMMON. 

AERODYNAMICS LOCAL COMMON TO GLOBAL. 
GLOBAL COMMON TO STABILITY LOCAL COMMON 


JCALC=1 
KCALC=2 
NERR = 0 


STEP 1. 

NSTRTR=ILOCAT (5,1) 

NMAXR=ILOCAT (5,2) 

NSTRTI=IL0CAT(5,3) 

CALL DATAIO (KCALC , NMAXR , NMAXI. , INOUT (NSTRTR) , INOUT(NSTRTI') , 
* RA5,IA5,RARAY, IARAY) 
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c 

c 


c 

c 

c 

c 


c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


STEP 2. 

NSTRTR=IL0CAT (3,1) 

NMAXR=IL0CAT(3,2) 

NSTRTI=IL0CAT(3,3) 

NMAXI=IL0CAT(3,4) 

CALL DATAIO (JCALC , NMAXR , NMAXI , INOUT (NSTRTR) , INOUT (NSTRTI) , 
* RA3 , IA3 , RARAY , I ARAY) 

STEP 3. 

CALL AER02 (ICALC , NERR , ICEO , KCPRNT ; IGPLT) 


STEP 4 . 

CALL DATAIO (KCALC, NMAXR, NMAXI, INOUT (NSTRTR) , INOUT (NSTRTI) , 
* RA3, IA3, RARAY, IARAY) 


STEP 5. 

IF (NERR . GT . 1) RETURN 
NSTRTR=ILOCAT (5,1) 

NMAXR=ILOCAT (5,2) 

NSTRTI=ILOCAT (5 , 3) 

NMAXI=IL0CAT(5,4) 

CALL DATAIO (JCALC, NMAXR, NMAXI , INOUT (NSTRTR) , INOUT (NSTRTI) , 
* RA5,IA5, RARAY, IARAY) 




THE FOLLOWING IS A LIST OF THE REFERENCES 
USED IN THIS SUBROUTINE : 

1) HOAK, D.E. et a I . ; "USAF Stability and Control DATCOM" ; 
Wright Patterson AFB, Ohio, 45433; Revised 1970 

2) Roskam, Jan; "Airplane Flight Dynamics and Automatic 
Controls, Part 1"; Published by the author, 519 Boulder 
Ave . Lawrence KA 66044; second printing 1979 

3) McCormick, B.W. ; "Aerodynamics, Aeronautics, and Flight 
Dynamics"; J. Wiley and Sons, 1979 

4) Perkins C.D. and Hage R.E.; "Airplane Performance, Stabil 
and Control"; J. Wiley and Sons, 1949 

5) Curry, Norman; "Aircraft Landing Gear Design, Principle 

and Practices"; ATAA Education Series; American Institute 
of Aeronautics and Astronaut i cs ; 1986 


i ty* 




******* 


RETURN 

END 
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THIS CONTAINS THE COMMON BLOCK STATEMENTS NEEDED TO RUN THE STABILITY 
AND CONTROL MODULE 'STBLCON . FOR ’ . THIS FILE IS CONNECTED WITH THE 
MODULE THROUGH THE 'INCLUDE' STATEMENT FOUND AT THE BEGINING OF EACH 
OF THE SUBROUTINES. 

COMMON/OVER/ I CALC , NERROR , MODNO , IGEO , KGPRNT , IGPLT , IF 

COMMON/STBCM/CLALFA , CLAFHT , DEDA , XQCBCN , XQCBWG , XQCBHT , XQCBVT , 

1 CBARCN , CBARWG , CBARHT , CBARVT , ROOTWG , TRWG , XLEWG , SPANWG , 

2 SWPWG , PODL , EN , BODL , XLEPOD (10) , WCAND , WFTOT , FUFRAC , WWING , WFS , 

3 WLG WHT , WVT , WETANK , WFEXT , WBODY , WAIRC , WAPU, WSC , WENG , SCN , SWG , SHT , 

4 SVT WCREW , WELT , WEP , WHDP , WINST , WNA , WPA , WCARGO , WAMMUN , WBOMB , WMISS , 

5 ARWG , DIHED , ZRTWC , ZRTHT , ARVT , SWPCN , TRCN , AMTO , 

6 ZRTCN, VOLB, BDMAX , FRN, FRAB, ALPHA , CL , THRUST , ARHT , TNT ( 1 2) , WFT ( 1 2 ) , 

7 WFTO WTOT , ARCN , CMWB , WGTO , WFEQ , WBB2 , WAF , WE , WPL , WPS , WTSUM , WFUEL , 

8 TWTO , MACH , ALT , SWETWC , SWETHT , SWETCN , SPANHT , SPANCN , 

9 STARTW(12) ,DLWMY, IAO 

COMMON /LOCALS/ETAHT , XFC AN , XFHT , XFNA , XFVT , XFWG , XFAPU , XFENG , 

1 XFELT , XFEP , XFINST , XFHDP , XFPA , XFPAYL , XFFFUS , XFCREW , WPAYL , XFFUR , 

2 XFLGFRT , WFLGFRT , XFLGR , XFFS , XFBOMB , XFSC , XFPI V , XFLIFTF , WPIV , WFUR , 

2 XEF1 , XFAMMUN , XFBB2 , XFMISS , XFCARGO, XACWG , ETACN , ODEDA , XAWB , CLACN , 

3 WZET , WZEF , WEFF , WEFW , XEF2 , WEFWT , WZENG , SUMMO , SIAAWT , XLE , LWB , LTTWB , 

3 IDBPRT,XCG(13) ,DCMDCL(13) ,SM(13) ,XCGAFT(14) , WLGFRT , WLGR , WFWG , 

4 WFP , WZB , WZM , WZ A , IBS , IMS , I AS , NPHASE , ICGPRT , WLFTF , WFFUS , LHT , LCN , 

5 C ELV , CBALHT , TZCHT , FNEUTPT ( 1 3 ) , XFBOD , XFFWG , XFEX , XLETIP , IETANK , 

6 ZMTOR , AWTOT , ZMU , XT1 , XT2 , ZFMG , ZMG , ZFT1 , ZT1 , ZFT2 , ZT2 , ZFD , ZD , Tl , 

7 TSPLIT , GAM ATI , GAMAT2 , IHT, ICN, ELVDEF , VROT, ALPHACN , TAU, T2 , 

8 SHTl , SCN1 , ELVDMAX , IWG , XCGFRWD (13) , QTOR , CGBAR (13) , ALFTO , CLTO , 

8 CDTO , CMWBTO , ALFMAX , CDMAX , CMWBMAX , ALFO , CLO , CDO , CMWBO , XLECN , 

9 Tl RAD , T2RAD , XLECG , XLET1 , XLET2 , CLMAX , CLAHTGRD , CLACNGRD , I VECT , 

1 SPANCN1 ,ZRTCN1 ,ARCN1 ,TRCN1 ,ALFHT r ALFCN,DCMTMUL,TWZTO,IFLAG 

COMMON/XCGS/XBOD , XWING , XCAN , XHT.XVT, XNA , XPIV , XAIRC , XAPU , 

1 XELT , XEP , XINST , XHDP , XSC , XPA , XENG (10) , XFS , XLIFTF , XFUR , 

2 XPAYL , XCREW , XCARGO , XAMMUN , XBB2 , XBOMB , XMISS , XETANK , 

3 XFUSF , XWGFUEL , XLGFRT , XLGR , LGFLAG 

REAL XACWG , DCMDCLF , DCMDCLT , DCMDCLW , DCMCLCN , KF , CSC , TC , MBOD , 

1 MCAN , MWG , MFS , MLC , MHT , MVT , MET ANK , MFEXT , MAIRC , MAPU , MSC , MENG (10) , 

1 LWB , LH1 , LH2 , LH3 , LH4 , LH5 , LH6 , LH7 , LH8 , LH9 , LC1 , LC2 , LC3 , LC4 , LC5 , 

2 LC6 , LC7 , LC8 , LHT , LCN , ZMTOR , ZMU , LTTWB , IWG , IHT , ICN , LEFF , MACH , 

2 MCREW , MELT , MEP , MHDP , MINST , MNA , MPA , MCARGO , MAMMUN , M8PMB , MMTSS , 

3 MPAYL ,MPIV,MLIFTF ,Ki 1 ,K!2,Ki?,Ki4 ,K21 ,K22.K23,K24 
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1 

0 115 

-61 

-62 

209 

100430 

100512 

100513 

211 

182 

529 

100069 

-307 

-308 

620 

619 

-610 

193 

-419 

421 

187 

568 

672 

673 

682 

683 

2 

0 392 

100013 

-60 

-241 

-242 

532 

-535 

100247 

100248 

100229 

100230 

100239 

100240 

100481 

100482 

330 

331 

451 

452 

271 

272 

312 

313 

322 

323 

443 

444 

263 

264 

-85 

-311 

100593 

100594 

68 

497 

289 

290 

299 

300 

401 

402 

217 

218 

227 

228 

561 

562 

491 

492 

625 

626 

159 

160 

120 

121 

130 

131 

22 

23 

151 

152 

100 

101 

20 

29 

569 

574 

347 

348 

360 

361 

379 

380 

501 

502 

170 

171 

180 

181 

100566 

694 


GEOMETRY REAL VARIABLE CATALOG 


-63 

-64 

100060 

100046 

100059 

100428 

100429 

100427 

100434 

100432 

100511 

-603 

-604 

-605 

-606 

196 

197 

282 

169 

100426 

73 

71 

72 

70 

370 

414 

412 

413 

411 

612 

663 

664 

634 

100066 

100607 

194 

195 

661 

184 

185 

418 

420 

417 

100422 

100388 

-567 

665 

666 

667 

668 

674 

675 

676 

677 

678 

684 

-570 

685 



TRAJECTORY 

REAL VARIABLE CATALOG 

-74 

-108 

-109 

100144 

100188 
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